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ABSTRACT 

This  thesis  investigates  the  accuracy  of  a  recently  proposed  passive  bearings-only 
Target  Motion  Analysis  (TMA)  procedure.  The  primar\'  method  of  analysis  is  to  com- 
pare computer  generated  positions  of  a  Target  that  is  moving  with  a  constant  course  and 
speed,  with  the  procedurally  derived  estimated  positions. 

A  computer  model  was  developed  which  simulated  several  possible  interactions  be- 
tween the  Target  and  Own  Ship.  Estimated  parameters  of  the  Target  track  were  com- 
puted using  the  procedure  under  analysis.  These  values  were  compared  to  the  "true" 
values  generated  by  the  simulation. 

An  analysis  of  the  TMA  procedure  as  originally  proposed  showed  that  it  failed  to 
accurately  estimate  the  target  track  parameters.  However  with  some  modifications,  the 
accuracy  improved  significantly  and  it  is  felt  that  the  procedure  can  accurately  estimate 
target  range  (but  not  necessarily  course  and  speed)  for  some  target  geometries. 
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I.     INTRODUCTION  TO  THE  PASSIVE  RANGING  PROBLEM 

A.     BRIEF  HISTORY 

The  importance  of  doing  Target  Motion  Analysis  (TMA)i  by  a  passive  bearings- 
only  tracking  technique  has  had  the  attention  of  many  tacticians  over  the  last  45  years, 
mainly  because  the  submarine's  objective  is  to  remain  covert  (undetected)  during  the 
tracking  maneuver.  This  guarantees  the  advantage  of  surprise  if  an  attack  is  to  be  per- 
formed. 

In  the  early  part  of  WWII,  TMA  by  passive  bearings-only  methods  began  to  be 
used,  and  the  following  drawn  from  Naval  Tactical  Decision  Aids  by  D.  H.  Wagner,  1989, 
provides  a  brief  histor\"  of  its  development: 

During  WWII  LT.  F.  C.  Linch  developed  the  Linch  Plot,  which  uses  a  pivotal 
relationship  among  bearings,  bearing  rate,  and  Target  relative  motion.  Based  on 
that,  he  developed  a  graphical  method  that  was  used  until  the  late  1960's.  It  was 
abandoned  because  it  was  useful  only  for  the  short  ranges  contemplated  by  Linch 
at  that  time. 

During  the  1950's,  various  human  plot  methods  and  nomographs  were  used  for 
TMA.  One  of  the  oldest  of  these  was  the  Strip  Plot,  later  called  Geographical  Plot. 
This  method  has  often  been  used  to  provide  upper  and  lower  bounds  on  Target 
range,  based  on  bounds  on  speed.  An  additional  much  used  nomographic  device 
has  been  the  Bearing  Rate  Slide  Rule.  An  early  TMA  computer  was  the  Position 
Keeper,  which  was  a  carryover  from  WWII,  and  used  as  an  aid  to  approach  and 
attack  surface  ships. 

In  1953.  F.  N.  Spiess  published  a  graphical  procedure  requiring  four  bearings  and 
a  maneuver  by  Own  Ship.  This  method  is  known  as  Spiess  Plot.  It  remains  in  use 
in  surface  ship  TMA. 

The  1954  command  thesis  of  LT.  J.  F.  Fagan  derived  a  four-bearing  TMA  sol- 
ution, which  was  a  transcendental  system  of  three  equations  with  three  unkno\^■ns. 
To  reduce  to  computability  by  slide  rule,  he  assumed  that  Own  Ship  motion  during 
the  first  three  bearings  was  approximately  zero  which  was  probably  satisfactory  for 
diesel  operations. 

In  1957,  LT.  J.  J.  Ekelund  devised  one  of  the  most  famous  TMA  methods,  now 

called  Lkelund  Ranging. 

Between  1967  and  1969,  D.  C.  Bossard  with  a  group  of  USX  ofTicers,  made  an 
analysis  of  the  fundamentals  of  bearings-only  ranging  and  the  elTect  of  Own  Ship 
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maneuvers  on  ranging  accuracy.  Through  this  work  and  trials  at  sea,  they  observed 
that  ranging  errors  could  be  eliminated  by  judicious  choice  of  the  time  for  which  the 
range  was  estimated.  Also,  by  selecting  some  particular  maneuvers,  this  best  time 
could  be  controlled  to  be  in  the  past,  present,  or  future;  then  they  applied  this 
method  to  various  forms  of  bearings-only  ranging,  especially  to  Ekelund  and  Spiess 
methods.   The  general  name  for  these  procedures  is  Time  Correction. 

In  August  196S,  a  series  of  exercises  was  conducted  aboard  the  USS  Pargo  to 
compare  the  different  methods  in  use  at  that  time.  The  Passive  Ranging  Manual 
evolved  in  three  volumes  from  the  report  of  the  analysis  of  these  exercises  and  was 
published  later  as  an  NWP.   This  is  a  basic  reference  today  for  TMA, 

In  the  late  1960's  the  CHURN  method  appeared,  which  uses  a  large  number  of 
bearing  observations  and  regression  methods  to  do  TMA.  This  method  was  devel- 
oped by  General  Dynamics'  Electric  Boat  and  Librascope. 

In  1969,  H.  \V.  Headle,  J.  Di  Russo,  and  E.  Messere  developed  the  Manual 
Adaptive  TMA  Evaluator  (MATE).  This  method  begins  with  an  input  estimate  of 
Target  course,  speed  and  range,  and  from  there  estimated  bearings  are  computed 
and  compared  against  observed  bearings.  DifTerences  between  these  bearing  values 
are  plotted  against  time.  If  these  differences  are  distributed  roughly  symmetrically 
as  white  noise  about  zero,  while  Own  Ship  changes  course  or  speed,  then  the  trial 
solution  is  a  good  one;  otherwise  the  MATE  operator  adjusts  the  input  Target  pa- 
rameters. 

In  the  early  1970's,  the  application  of  Kalman  filtering  appeared  for  the  first  time 
in  TMA.  It  was  a  done  by  J.  S.  Davis  of  XUSC,  based  on  the  Ph.  D.  thesis  of  D. 
J.  Murphy  at  Northeastern  University.    Later  versions  were  developed  by  IBM. 

In  the  mid  1970's,  a  new  TMA  method  known  as  FLIT  was  devised  by  ENS.  L. 
Anderson  which  has  seen  considerable  elTective  operational  use.  Its  methodology 
remains  classified. 

In  the  mid  and  later  1970's,  ranging  and  tracking  at  long  range  on  a  sphere  were 
developed  by  D.  C.  Bossard,  J.  B.  Behrla.  and  L.  K.  Graves  of  Daniel  H.  Wagner, 
Associates.    This  was  motivated  bv  over-the-horizon  targeting  needs. 
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In  the  late  1970's  and  during  the  1980's.  further  developments  were  made  by  using 
Kalman  filtering  and  also  stochastic  difierential  equations.  The  Maneuvering  Target 
Statistical  Tracker  (MTST)  is  a  method  devised  from  this  basis  by  W.  H.  Barker  of 
Daniel  H.Wagner,  Associates. 

From  1985  to  19cS7  the  Generical  Statistical  Tracker  (GST)  was  developed  at 
Daniel  H.  Wagner,  Associates.    This  is  a  PC  version  of  MTST. 

B.     PURPOSE  OF  THE  THESIS 

In  late   198S,  LCDR  P.   K.  Peppe  of  the  USS  La  Jolla  (SSN  701)  proposed  a 
bearmg-only  Target  Motion  Analysis  (TMA)  procedure  (here  called  the  Bearings  Ex- 


trapolation  Procedure)  for  use  in  submarines  [Ref.  1].    It  is  the  purpose  of  this  thesis  to 
test  the  accuracy  of  this  procedure  against  a  nonmaneuvering  Target. 

C.     PROBLEM  DESCRIPTION 

A  submarine  (Own  Ship)  is  tn,'ing  to  determine  the  range,  course,  and  speed  of  a 
ship  (Target)  that  is  moving  with  a  fixed  course  and  speed.  The  data  collected  by  Own 
Ship  consists  of  passive  bearings.  Since  the  bearing  data  is  obtained  by  a  single  tracker, 
a  unique  tracking  solution  may  only  be  obtained  by  a  maneuvering  Own  Ship  [Ref  2]. 
A  TMA  method  which  uses  a  two-leg  maneuver  can  fulfill  this  requirement. 

The  two-leg  TMA  procedure  used  in  this  study  is  depicted  in  Figure  1.  Own  Ship 
begins  by  receiving  bearings  to  the  Target  each  20  seconds  for  4  minutes.  In  Figure  I, 
t,  is  the  time  at  which  bearing  /  is  received.  After  time  t^^,  the  bearings  stop  and  Own 
Ship  begins  a  course  change.  At  time  /3c,  Own  Ship  is  on  the  second  leg  and  bearings 
resume.  Bearings  continue  on  the  second  leg,  the  last  one  being  received  at  time  /4,. 
Estimates  of  Target  course,  speed  and  range  will  be  derived  from  these  24  observed 
bearings. 

The  central  idea  of  the  Bearings  Extrapolation  Procedure  is  that  given  three  bearings 
observed  at  three  times  on  a  single  leg.  all  future  and  past  bearings  can  theoretically  be 
determined.  (This  idea  was  known  by  Spiess  [Ref  3]  and  may  predate  him.)  So,  refer- 
ring to  Figure  1,  the  estimated  Target  position  at  time  t^  can  be  obtained  by  crossing  the 
true  bearing  (beginning  at  point  /(,)  with  the  extrapolated  bearing  from  leg  two  (begin- 
ning at  point  /fl'). 

The  same  as  above  can  be  done  for  all  time  points  that  have  bearings  between  i^  and 
^4,.  Then  regressing  linearly  on  these  estimated  Target  positions,  the  Target  track  pa- 
rameters can  be  determined. 

Although  this  procedure  should  give  exact  answers,  in  practice  large  errors  can  be 
generated.  The  reason  for  these  inaccuracies  is  that  bearings  to  the  Target  cannot  be 
measured  exactly,  so  the  extrapolated  bearings  will  also  be  in  error. 
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Figure  I.      Graphical  Representation  of  the  Bearings  Extrapolation  Procedure. 


Exactly  how  bearings  that  contain  errors  are  extrapolated,  with  as  much  accuracy 
as  possible,  forward  and  backward  in  time  is  part  of  what  this  thesis  examines.  The  most 
straightforward  method  (and  the  first  to  be  analyzed)  is  simply  to  use  linear  regression 
on  the  {B,  ,  t)  data,  where  B,  is  the  bearing  at  time  t.  Alternate  methods  include  re- 
gressing the  {B, ,  t)  data  by  using  quadratic  or  arctangent  models.  Each  of  these  proce- 
dures were  tried  and  accuracies  assessed. 


The  Spiess  Graphical  Method,  Ekelund  Ranging,  and  the  CHURN  Method  are 
closely  related  to  the  procedure  proposed  by  Peppe,  and  for  that  reason  they  are  further 
explained  in  Appendix  A. 

D.  ASSUMPTIONS  FOR  THE  SIMULATION 

The  computer  simulation  used  in  this  study  uses  the  following  assumptions: 

1.  The  Target  moves  at  a  constant  course  and  speed  throughout  the  tracking  maneu- 
ver. 

2.  Differences  in  depth  for  the  Target  and  Own  Ship  are  disregarded. 

3.  Errors  in  position  for  Own  Ship  are  disregarded. 

4.  Errors  in  bearings  to  the  Target  are  independent  and  identically  distributed  normal 
random  variables  with  mean  zero  and  standard  deviation  not  larger  than  1.5°.  In 
addition,  this  standard  deviation  is  constant  for  each  simulation  experiment. 

5.  Bearings  to  the  Target  are  obtained  each  20  seconds. 

6.  Initial  course,  speed  and  range  for  the  Target  and  initial  course  and  speed  for  Own 
Ship  are  as  described  in  Chapter  II. 

E.  TERMINOLOGY 

1.  Lead-leg.  An  Own  Ship  course  such  that  Own  Ship  speed  across  the  line  of  sound 
is  in  the  same  direction  as  the  Target  speed  across  the  line  of  sound  at  the  start  of 
the  leg. 

2.  Lag-leg.  An  Own  Ship  course  such  that  Own  Ship  speed  across  the  line  of  sound 
is  in  the  opposite  direction  to  the  Target  speed  across  the  hne  of  sound  at  the  start 
of  tlie  leg. 

3.  Leg  length.  The  distance  along  either  leg  where  bearings  are  received  (does  not 
include  Own  Ship  turn). 

4.  Raw  bearing.   An  unsmoothed  bearing  to  the  Target. 

5.  Faired  bearing.   A  bearing  obtained  after  smoothing  of  the  raw  bearings. 

6.  Extrapolated  positions.  Positions  of  Own  Ship  obtained  by  extending  the  first  leg 
forward  in  time  or  the  second  leg  backward  in  time. 

7.  Extrapolated  bearings.  Bearings  obtained  by  extrapolating  actual  bearings  forward 
in  time  for  the  first  leg  or  backward  in  time  for  the  second  leg. 

8.  Lead-leg  estimation.  The  estimated  target  track  obtained  by  linear  regression  of 
estimated  Target  positions.  These  estimated  Target  positions  are  obtained  from  the 
intersection  of  faired  bearings  from  the  lead-leg  and  extrapolated  bearings  from  the 
lag-leg. 

9.  Lag-leg  estimation.  The  estimated  target  track  obtained  by  linear  regression  of 
estimated  Target  positions.  These  estimated  Target  positions  are  obtained  from  the 
intersection  of  faired  bearings  from  the  lag-leg  and  extrapolated  bearings  from  the 
lead-leg. 


10.  Lead  angle.   Angle  between  the  line  of  sound  and  the  course  of  Own  Ship  when  it 
is  leading  the  Target. 

1 1.  Lag  angle.   The  angle  between  the  line  of  sound  and  the  course  of  Own  Ship  when 
it  is  lagging  the  Target. 

12.  Angle  on  the  bow.    The  angle  between  the  line  of  sound  and  the  course  of  the 
Target. 


II.     TESTING  THE  BEARING  EXTRAPOLATION  PROCEDURE 

A.  INITIAL  CONDITIONS 

In  order  to  perform  a  thorough  analysis  of  the  Bearing  Extrapolation  Procedure, 
5,184  separate  simulation  experiments  were  conducted.  Each  experiment  used  one  of  the 
possible  combinations  of  the  following  parameters: 

1.  The  initial  Target  positions  were  all  due  north  from  Own  Ship  (i.  e.,  000°)  at  one 
of  three  different  initial  distances:  large  (=  60,000  yds),  medium  (=  30,000  yds), 
and  short  (=   10,000  yds). 

2.  Six  courses  for  the  Target  were  used:   030°,  060°,  090°,  120°,  150°,  and  175°. 

3.  Three  speeds  for  the  Target:   25,  15,  and  5  knots. 

4.  Four  courses  for  Own  Ship  in  the  lead-leg:   060°,  070°,  080°,  and  090°. 

5.  Four  courses  for  Own  Ship  in  the  lag-leg:    270°,  280°,  290°,  and  300°. 

6.  Three  speeds  for  Own  Ship:   5,  10,  and  15  knots. 

7.  Two  maneuvers  by  Own  Ship:   lead-lag  and  lag-lead. 

8.  Four  bearing  error  standard  deviations:  0°,  0.5°,  1°,  and  1.5°. 

It  is  important  to  mention  that  because  of  problem  symmetn.'  and  the  fact  that  both 
lead-lag  and  lag-lead  maneuvers  are  simulated,  there  was  no  need  for  simulating  courses 
for  the  Target  from  180°  to  36(r. 

The  randomness  in  the  simulation  was  introduced  only  through  errors  in  the  re- 
ceived bearings. 

B.  DESCRIPTION  OF  THE  SIMULATION 

The  simulation  was  written  in  Fortran  77  [Ref  4]  and  executed  on  an  IBM  3033 
mainframe.  Appendix  B  has  a  complete  listing  of  all  variables  of  the  simulation  pro- 
gram.  Appendix  C  is  a  listing  of  the  simulation  program  code. 

For  each  combination  of  initial  conditions,  the  simulation  was  repeated  100  times 
with  different  seeds  for  pseudorandom-generation  of  bearing  error.  For  each  of  the  re- 
petitions, the  following  steps  were  conducted: 

1.  Target  positions  were  simulated  every  20  seconds  for  the  entire  maneuver. 

2.  For  the  first  leg  the  following  data  was  generated: 

a.  Own  Ship  position  every  20  seconds. 

b.  Bearing  to  the  Target  every  20  seconds  (normal  errors  are  added). 


c.   Own  Ship  track  extended  beyond  the  turn. 

3.  For  the  first  leg  the  following  data  was  then  computed: 

a.  CoefTicients  for  linear  and  quadratic  fits  for  bearings  versus  time. 

b.  Faired  bearings  and  extrapolated  bearings  using  the  results  from  linear  or 
quadratic  fit. 

4.  Own  Ship  course  was  determined  for  the  second  leg;  based  on  the  required  second 
leg,  lead  or  lag  angle,  and  the  line  of  sound  specified  by  the  extrapolated  bearing 
at  the  middle  of  the  maneuver  (/^o). 

5.  Simulations  and  computations  for  the  second  leg  w^ere  completed  in  the  same 
fashion. 

6.  Estimated  positions  of  the  Target  in  two  coordinates  (X  and  Y)  were  computed  by 
determining  the  point  of  intersection  of  faired  bearings  with  their  corresponding 
extrapolated  bearings. 

7.  Linear  and  orthogonal  regression  is  conducted  on  the  estimated  Target  positions 
(in  X  versus  ^'  coordinates),  to  determine  the  Target  course. 

8.  Range  to  the  Target  at  the  end  of  the  maneuver  was  obtained  by  computing  the 
distance  between  Own  Ship's  actual  position  at  the  end  of  the  second  leg  and  the 
estimated  Target  position. 

9.  Finally.  Target  speed  was  obtained  by  computing  the  distance  betw'een  the  esti- 
mated Target  position  at  time  ^o  and  the  estimated  Target  position  at  the  end  of  the 
second  leg  (/j,),  on  the  fitted  Target  track,  and  then  dividing  this  distance  by 

'41    -    /O- 

The  simulation  of  the  bearing  errors  was  done  using  the  Linear  Congruential 
Method  for  uniform  pseudorandom  numbers,  combined  with  the  Box  and  Muller 
Method  to  produce  normal  variates  [Refs.  5,6].  These  procedures  were  selected  to  allow 
the  simulation  to  be  performed  on  a  personal  computer. 

C.     TESTLNG  THE  BEARING  EXTRAPOLATION  PROCEDURE  UNDER  IDEAL 
CONDITIONS 

Flere  bearings  to  the  Target  were  considered  without  error.  Also  the  bearing  rates 
immediately  before  the  turn  (used  to  extrapolate  bearings  forward  in  time)  and  imme- 
diately after  the  turn  (used  to  extrapolate  bearings  backward  in  time)  were  computed 
exactly  with  no  errors  introduced.  The  only  source  of  error  in  this  case  was  the  linear 
bearing  extrapolation.  This  case  is  intended  to  be  an  optimistic  bound  on  the  accuracy 
of  the  bearing  extrapolation  procedure  when  linear  regression  is  used. 

The  measure  of  efiectiNcness  (MOE)  for  the  accuracy  of  this  procedure  was  the 
percentage  of  times  that  the  estimated  error  was  within  10"o  for  range.  20°  for  course, 
and  10",,  for  speed  from  the  actual  parameters  of  the  Target  track.    These  results  are 


presented  in  Table  1  and  Table  2  for  lead-lag  and  lag-lead  maneuvers  respectively. 
Each  number  in  these  Tables  resulted  from  the  examination  of  18  simulation  exper- 
iments (six  Target  courses  and  three  Target  speeds).  For  example,  when  the  initial  dis- 
tance was  large,  lead  and  lag  angles  70°,  and  Own  Ship  speed  15  knts,  then  five  of  the 
18  simulations  experiments  (or  27.8%)  resulted  in  the  calculated  Target  speed  being 
within  10%  of  the  actual  Target  speed.  Note  that  without  bearing  error,  the  simulation 
experiments  are  deterministic. 


Table  I.     ACCURACY  IN  IDEAL  CONDITIONS  FOR  BEARINGS  EXTRAPO- 
LATION PROCEDURE  (LEAD-LAG  MANELIATR) 


Own  Ship 
speed 
(knts) 

Lead  and 
lag  angle 

C) 

Initial  dis- 
tance 

Percentage  of  times  the  absolute  error  is 

smaller  than  10°  o  of  actual  for  range  and 

speed  or  20°  of  actual  for  course. 

Range 
(10%) 

Course 
(20°) 

Speed 
(10%) 

15 

70 

Large 

100 

50.0 

27.8 

Medium 

100 

50.0 

27.8 

Short 

38.9 

50.0 

22.2 

SO 

Large 

100 

50.0 

33.3 

Medium 

100 

55.6 

33.3 

Short 

50.0 

50.0 

27.8 

90 

Large 

100 

55.6 

27.8 

Medium 

100 

55.6 

27.8 

Short 

50.0 

66.7 

33.3 

10 

70 

Large 

100 

44.4 

44.4 

Medium 

94.4 

38.9 

-1  -s    -> 

Short 

55.6 

38.9 

22.2 

80 

Large 

100 

44.4 

33.3 

Medium 

94.4 

44.4 

33.3 

Short 

55.6 

44.4 

27.8 

90 

Large 

100 

55.6 

33.3 

Medium 

94.4 

50.0 

Short 

61.1 

50.0 

27.8 

5 

70 

Large 

83.3 

44.4 

22.2 

Medium 

72.2 

44.4 

27.8 

Short 

44.4 

44.4 

27.8 

80 

Large 

88.9 

33.3 

27.8 

Medium 

77.8 

33.3 

27.8 

Short 

44.4 

44.4 

->-»  -> 

90 

Large 

8S.9 

2~.8 

Medium 

77.8 

33.3 

">")    T 

Short 

44.4 

33.3 

16.7 
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Table  2.     ACCURACY  IN  IDEAL  CONDITIONS  FOR  BEARINGS  EXTRAPO- 
LATION PROCEDURE  (LAG-LEAD  MANELFVER) 


Own  Ship 
speed 
(knts) 

Lead  and 
lae  ancle 

Initial  dis- 
tance 

Percentage  of  times  the  absolute  error  is 

smaller  than  lO'^/o  of  actual  for  range  and 

speed  or  20°  of  actual  for  course. 

Ranee 

(I0"o) 

Course 
(20°) 

Speed 

(10°;) 

15 

70 

Large 

88.9 

77.8 

27.8 

Medium 

66.7 

77.8 

27,8 

Short 

50.0 

50.0 

33.3 

80 

Large 

8S.Q 

94.4 

38.9 

Medium 

72.2 

88.9 

50.0 

Short 

50.0 

33.3 

33.3 

90 

Large 

94.4 

94.4 

44.4 

Medium 

72.2 

88.9 

44.4 

Short 

38.9 

27.8 

27.8 

10 

70 

Large 

83.3 

88.9 

22.2 

Medium 

66.7 

94.4 

")■>  0 

Short 

44.4 

50.0 

27.8 

80 

Large 

83.3 

94.4 

Medium 

61.1 

88. 8 

33.3 

Short 

50.0 

44.4 

1  -I  -> 

90 

Large 

77.8 

94.4 

38.9 

Medium 

61.1 

88.9 

38. 9 

Short 

50.0 

38.9 

16.7 

5 

70 

Large 

72.2 

77.8 

3 -5.  J 

Medium 

50.0 

77.8 

27.8 

Short 

3S.^) 

66.7 

27.8 

80 

Large 

■7-)  -) 

77.8 

33.3 

Medium 

50,0 

77.8 

27.8 

Short 

44,4 

66.7 

27,8 

90 

Large 

72.2 

77.8 

38.9 

Medium 

50. n 

77.8 

Short 

44.4 

66.7 

->■>  -) 
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As  can  be  seen  in  Table  1  and  Table  2,  the  following  results  were  obtained  under 
ideal  conditions: 

1.  Range  estimation  is  good,  but  only  for  medium  or  large  initial  distances. 

2.  Accuracy  improves  with  increases  in  Own  Ship  speed. 

3.  Course  and  speed  estimation  is  less  accurate  than  range  estimation. 

4.  Accuracy  improves  with  increases  in  lead  and  lag  angle. 

5.  Lead-lag  maneuvers  give  better  range  estimates  than  do  lag-lead,  but  not  neces- 
sarily better  course  and  speed  estimates. 

It  is  felt  that  the  procedure  performed  poorly  for  short  initial  distances,  in  part  be- 
cause the  bearing  versus  time  curve  is  more  nonlinear  in  this  case. 

D.     TESTING  WITH  BEARING  ERRORS 

For  these  experiments,  the  program  was  run  with  bearing  errors  added  to  the  simu- 
lated bearings.  Also  the  bearing  rate  for  each  leg  was  computed  using  a  linear  least 
squares  technique. 

Only  lead  or  lag  angles  of  90°  were  considered,  since  it  was  determined  in  the  pre- 
vious step  that  this  is  the  best  course  for  Own  Ship  to  follow.  Accuracy  in  course  and 
speed  were  not  considered  since  the  errors  for  these  parameters  were  too  large,  even  in 
ideal  conditions.  An  Own  Ship  speed  of  5  knots  was  not  considered  for  the  same  reason. 

The  MOE  for  these  experiments  was  the  percentage  of  times  (runs)  that  the  esti- 
mated range  at  time  /ji  was  within  5°^,  10%.  and  20%  of  the  actual  range.  Each  number 
in  Table  3  and  Table  4  results  from  1,800  trials  (six  Target  courses,  three  Target  speeds 
and  100  replications  per  (course,  speed)  combination). 
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Table  3.     ACCURACY  IN   REAL  CONDITIONS   FOR  BEARING  EXTRAPO- 
LATION PROCEDURE  (LEAD-LAG  MANEirv  ER) 


Own  Ship 
speed 
(knts) 

St.  Dev. 

bearing  er- 
ror 

Initial  dis- 
tance 

Percentage  of  times  the  absolute  error  is 
smaller  than  5"/o,  10%,  and  20%  of  ac- 
tual for  range. 

^       0 

10% 

20"o 

15 

1.5° 

large 

6.2 

11.8 

25.5 

medium 

12.3 

24.4 

45.4 

short 

15.9 

30.6 

53.7 

1.0° 

large 

8.9 

18.4 

36.6 

medium 

15.8 

37.8 

60.3 

short 

18.7 

38.3 

60.4 

0.5° 

large 

16.9 

33.7 

61.9 

medium 

26.5 

52.5 

82.2 

short 

19.4 

39.2 

62.8 

10 

1.5° 

large 

4.2 

8.8 

18.2 

medium 

8.4 

17.0 

33.9 

short 

12.6 

24.4 

51.2 

1.0° 

large 

5.8 

12.3 

26.4 

medium 

10.7 

23.5 

45.3 

short 

15.1 

30.2 

59." 

0.5° 

large 

lo.S 

22.9 

45.  S 

medium 

19.0 

38.9 

68.7 

short 

20.3 

38.8 

69.4 
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Table  4.     ACCURACY  IN   REAL   CONDITIONS   FOR  BEARING  EXTRAPO- 
LATION PROCEDURE  (LAG-LEAD  MANEUV  ER) 


Own  Ship 
speed 

(knts) 

St.  Dev. 

bearing  er- 
ror 

Initial  dis- 
tance 

Percentage  of  times  the  absolute  error  is 
smaller  than  5%.  10%,  and  20%  of  ac- 
tual for  range. 

5% 

10% 

20% 

15 

1.5° 

large 

4.7 

11.3 

23.2 

medium 

10.3 

19.7 

36.5 

short 

6.7 

12.3 

23.9 

1.0° 

large 

7.11 

15.9 

34.3 

medium 

11.7 

24.2 

46.6 

short 

6.7 

12.2 

22.2 

0.5° 

large 

16.2 

31.8 

56.1 

niedium 

17.1 

30.9 

54.9 

short 

5.3 

11.6 

21.4 

10 

1.5° 

large 

4.0 

7.4 

18.1 

medium 

7.3 

13.7 

26.9 

short 

5.7 

11.7 

21.9 

1.0° 

large 

5.4 

10.8 

23.4 

medium 

9.(1 

17.3 

35.9 

short 

5.0 

10.8 

23.2 

0.5° 

large 

11.(1 

7")  ") 

42.9 

medium 

12.4 

25.1 

46.2 

short 

3.8 

9.7 

23.6 

From  Table  3  and  Table  4,  the  following  results  are  obtained: 

1 .  The  higher  the  Own  Ship  speed,  the  better  the  accuracy. 

2.  The  best  accuracy  is  obtained  for  medium  initial  distances. 

3.  The  smaller  the  bearing  error,  the  better  the  accuracy. 

4.  A  lead-lag  maneuver  gives  better  accuracy  than  lag-lead  maneuver. 

It  is  noted  that  orthogonal  regression  was  also  tried  on  the  estimated  Target  posi- 
tions, but  with  results  similar  to  those  obtained  with  linear  regression.  Subsequent  ex- 
periments used  only  linear  regression. 
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E.     ANALYSIS  OF  THE  RESULTS 

It  appears  that  even  in  the  absence  of  bearing  errors,  the  Bearing  Extrapolation 
Procedure  gives  poor  estimates  of  Target  course  and  speed.  And  when  bearing  errors 
are  introduced,  the  range  errors  also  become  unacceptable.  It  is  the  purpose  of  this 
section  to  determine  why  these  poor  results  were  obtained. 

L     Differences  Between  Lead-Leg  and  Lag-Leg-Estimation 

Figures  2  to  5  show  actual  and  estimated  Target  positions  for  four  encounter 
geometries  with  no  error  in  bearing.  So  the  only  errors  introduced  are  due  to  the  linear 
extrapolation  of  bearings.  It  is  clear  that  lag-leg  estimation  is  much  more  accurate  than 
the  lead-leg-estimation.  It  is  also  clear  that  any  attempt  to  regress  estimated  Target 
positions  over  both  legs  can  result  in  large  errors. 
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Figure  2.       Actual  (.)  and  Estimated  (*)  Target  Track  for  Lead- Lag  Maneuver. 
(Target  Course  =  060=) 


15 


ACTUAL  TARGET  7KACK 


\  •■' 


LAC  LEG  ESTIMATION 


f 


LEAD  LEG  ESTIMATION 


2000  40CX) 

TARGET  POSrriON  IN  Y  (y<j») 


sooo 


Figure  3.       Actual  (.)  and  Estimated  (*)  Target  Track  for  Lead-Lag  Maneuver. 
(Target  Course  =  150°) 
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Figure  4.       Actual  (.)  and  Estimated  (*)  Target  Track  for  Lag-Lead  Maneuver. 
(Target  Course  =  030°) 
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Figure  5.       Actual  (.)  and  Estimated  (*)  Target  Track  for  Lag-Lead  Maneuver. 
(Target  Course  =   120°) 

2.     Explanation  of  Why  Lag-Leg  Estimation  is  Better  than  Lead-Leg  Estimation 

To  help  explain  why  the  lag-leg  estimation  was  better  than  lead-leg  estimation, 
scatter  plots  of  bearings  versus  time  for  numerous  encounter  geometries  were  obtained. 
Figures  6,  7,  and  8  are  typical.  In  Figure  6,  a  lead-lag  encounter  is  represented.  It  was 
observed  that  during  the  lead-leg  (first  leg),  the  bearing  rate  was  generally  smaller  than 
that  of  the  lag-leg  (second  leg).  Also,  and  perhaps  more  importantly,  the  change  in 
bearing  rate  (bearing  acceleration)  was  also  smaller  in  the  lead-leg.  When  this  happens, 
it  is  possible  to  use  linear  regression  to  extrapolate  the  lead-leg  bearings  accurately  into 
the  future.  This  makes  the  lag-leg  estimation  more  accurate.  It  is  clear  from  Figure  6 
that  if  linear  regression  were  used  to  extrapolate  the  lag-leg  bearings  into  the  past,  sig- 
nificant errors  would  result,  and  these  errors  would  likely  cause  the  lead-leg  estimation 
to  be  less  accurate. 

In  the  same  fashion.  Figure  7  shows  a  plot  of  bearings  versus  time  for  a  typical 
lag-lead  encounter.  Here  also,  the  lead-leg  (second  leg)  has  a  smaller  bearing  rate  and 
bearing  acceleration.  So  linear  regression  can  be  used  to  extrapolate  accurately  the 
lead-lee  beannss  backward  m  time,  acain  makini:  the  lae-lee  estimation  most  accurate. 
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Unfortunately,  however,  there  are  encounter  geometries  where  neither  the  lead 
nor  the  lag-legs  show  small  bearing  acceleration.  For  example,  Figure  8  shows  a  case 
where  the  Target  is  at  high  speed,  short  initial  distance,  and  on  course  160°.  Here  both 
lead  and  lag-leg  bearings  are  nonlinear,  making  the  Bearing  Extrapolation  Procedure 
unusable.  Extensive  experimentation  indicates  that  when  all  the  following  conditions 
hold: 

Target  course  >  150  °, 

Target  speed  >  5  knots, 

Target  range   >  40,000  yds; 

then  the  lead-leg  bearings  (and  the  lag-leg  bearings)  are  not  significantly  linear  to  be 
accurately  extrapolated,  and  as  a  result,  lag-leg  estimation  (as  well  as  lead-leg  esti- 
mation) is  poor. 

Finally,  it  is  noted  that  for  very  long  initial  distances  (>  40,000  yds),  both  the 
lead  and  lag-legs  have  approximately  a  linear  change  in  the  bearings,  making  both  lead 
and  las-lee  estimations  reasonably  accurate. 
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Figure  6.      Bearings  in  Both  Legs  for  Lead-Lag  Maneuver  (Target  course  =  060°) 
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Figure  7.      Bearings  in  Both  Legs  for  Lag- Lead  Maneuver  (Target  Course  =   140°) 
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Figure  8.       Bearings  in  Both  Legs  for  Lead-lag  Maneuver  at  Sliort  Distance  and 
Closing  Target. 
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F.     CONCLUSIONS 

1.  Even  when  bearings  were  without  error,  the  Bearing  Extrapolation  Procedure  was 
unsuccessful  in  estimating  Target  course  and  speed.  However,  range  estimates  were 
reasonable. 

2.  When  bearing  errors  were  introduced,  range  estimates  also  became  inaccurate. 

3.  For  most  encounter  geometries,  lag-leg  estimation  gave  more  accurate  results  than 
did  lead-leg  estimation.  However,  for  some  important  cases  both  lag  and  lead-leg 
estimations  performed  poorly. 
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III.     POSSIBLE  IMPROVEMENTS,  AND  ANALYSIS  OF  THEIR 

ACCURACIES 

The  lack  of  accuracy  of  the  proposed  procedure,  even  under  ideal  conditions,  led  to 
attempts  to  improve  its  performance.  The  following  are  some  suggested  minor  changes 
and  two  suggested  major  changes  that  were  found  to  improve  the  accuracy  of  the  Bear- 
ing Extrapolation  Procedure. 

A.     MINOR  CHANGES 

1.  Increase  in  Leg  Length 

In  an  attempt  to  improve  performance,  several  different  leg  lengths  were  tried. 
It  was  discovered  that  leg  lengths  of  S  to  12  minutes  gave  the  best  results.  Shorter  legs 
gave  too  little  data  and  longer  legs  resulted  in  nonUnear  changes  in  lead-leg  bearings. 
Based  on  these  tests,  ten  minutes  was  selected  as  the  leg  length  for  subsequent  simu- 
lation experiments. 

2.  Use  of  a  Quadratic  Model  to  Fair  the  Lag-Leg 

Another  important  change  from  the  original  simulation  that  was  tried  was  to 
use  a  quadratic  model  to  fair  the  bearings  in  the  lag-leg.  This  produced  better  accuracy 
in  the  estimated  Target  track  parameters.  It  is  noted  that  current  manual  TMA  methods 
(in  particular,  the  Time-Bearing  Plot)  use  linear  fairing. 

3.  Increase  in  Turn  Time 

Contrary'  to  Ekelund  ranging,  if  an  instantaneous  turn  is  considered,  large  errors 
in  the  estimated  Target  track  parameters  are  generated.  This  is  because  some  actual  and 
corresponding  extrapolated  bearings  are  too  close  to  each  other  (i.  e.,  a  small  baseline). 
The  intersection  of  those  bearings  often  results  in  a  estimated  Target  position  far  away 
from  the  true  position  of  the  Target. 

Based  on  these  considerations,  it  was  determined  that  the  time  required  for  Own 
Ship  to  turn  should  be  increased.  After  many  trials,  it  was  found  that  the  best  turn  du- 
ration is  ten  minutes,  if  leg  length  is  also  ten  minutes,  For  longer  times,  the  bearing  rate 
in  the  lead-leg  does  not  remain  low  and  constant,  reducing  the  accuracy  of  the  lead-leg 
bearing  extrapolation. 
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B.     ELIMINATION  OF  THE  LEAD-LEG-ESTIMATION 

As  shown  in  Chapter  II,  when  linear  extrapolation  of  the  bearings  from  the  lag-leg 
is  used,  increased  errors  are  generated  in  the  estimation  of  the  Target  track  parameters. 
This  occurs  because  the  lag-leg  generates  a  high  and  variable  bearing  rate  compared  to 
the  lead-leg.  A  simple  solution  to  the  problem  is  to  drop  the  estimated  Target  positions 
generated  by  the  extrapolation  of  the  lag-leg  bearings. 

One  problem  with  this  idea  is  that  a  lead-lag  maneuver  will  then  produce  more  ac- 
curate estimates  of  Target  position  at  the  end  of  the  maneuver  than  will  a  lag-lead  ma- 
neuver. This  occurs  since  in  the  lead-lag  case,  the  estimated  Target  positions  are 
obtained  for  times  during  the  second  (i.e.,  the  lag)  leg.  And  in  the  lag-lead  case,  the  es- 
timated Target  positions  are  for  times  in  the  first  (again,  the  lag)  leg.  So  to  use  a  lag-lead 
maneuver  to  estimate  Target  positions  at  the  end  of  the  maneuver  requires  that  these 
positions  be  extrapolated  in  time  using  derived  estimates  of  Target  course  and  speed. 
This  leads  to  accumulated  errors  in  the  final  Target  position. 

To  test  for  the  possible  improvements  that  might  result  when  using  only  lag-leg  es- 
timation, simulation  experiments  were  conducted  with  the  following  parameters: 

1.  Speed  of  Own  Ship:   5,  10,  and  15  knots. 

2.  Speed  of  the  Target:   5,  15.  and  25  knots. 

3.  Initial  range:    large  (  =  60.000),  medium  ( =  30,000),  and  short  ( =  10,000)  yds. 

4.  Only  the  lead-lag  maneuver  is  considered. 

5.  Only  the  initial  lead  angle  of  90'  is  cosidered. 

6.  Leg  length:   4  minutes. 

7.  Time  between  legs:    6  minutes. 

8.  Standard  deviation  bearing  error:   0°,  0.5°.  1.0°,  and  1.5°. 

9.  Bearings  for  the  first  leg  were  faired  and  extrapolated  linearly. 

10.  Bearings  in  the  second  leg  were  faired  linearly. 

It  is  important  to  note  that  other  values  for  lead  angle  were  not  considered  because, 
as  discussed  in  Chapter  11,  it  was  found  that  the  best  results  were  obtained  for  a  90°  lead 
angle. 

In  order  to  show  how  the  minor  changes  improve  the  performance  of  this  major 
suggested  change  (i.  e..  elimination  of  the  lead-leg  estimation),  the  simulation  was  run 
for  three  different  cases: 
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1.  Without  Minor  Changes 

The  simulation  program  was  run  100  times  for  each  combination  of  the  above 
parameter  values.  It  was  determined  that  for  the  ideal  case  (0°  bearing  errors)  the  ac- 
curacy of  the  predicted  Target  parameters  was  adequate,  but  with  bearing  errors  intro- 
duced, large  errors  in  all  three  Target  parameters  were  produced.  So  it  was  concluded 
that  simply  eliminating  the  lead-leg  estimation  was  not  a  sufficient  improvement. 

2.  With  Increased  Leg  Length  and  Lag-Leg  Faired  with  a  Quadratic  Model 

Here  leg  length  was  increased  to  ten  minutes  and  the  lag-leg  was  faired  with  a 
quadratic  model.  The  same  simulation  experiments  were  run  as  before.  Only  the  results 
for  lead-lag  maneuvers  are  presented  in  Table  5.  Lag-lead  maneuvers  gave  uniformly 
worse  results  for  reasons  discussed  above.  Additionally,  and  more  surprisingly,  the  ac- 
curacy for  the  lag-lead  maneuver  was  not  as  good  as  the  lead-lag,  even  when  estimating 
the  Target  parameters  at  the  end  of  the  first  leg.  This  was  a  confirmation  that,  in  gen- 
eral, the  lead-lag  maneuver  gives  better  results  than  do  lag-lead  [Ref  7:  pp.  3]. 

Results  for  estimated  Target  course  and  speed  were  very  precise  only  when  the 
bearing  error  was  zero  and  when  Own  Ship  speed  was  greater  than  ten  knots.  When 
errors  in  bearings  were  included,  the  accuracy  in  estimated  Target  course  and  speed  was 
tremendously  reduced.  For  example,  for  large  initial  distance  and  15  knots  Own  Ship 
speed  the  following  was  obtained: 

1.  With  no  bearing  error,  the  error  in  estimated  speed  was  always  less  than  W^/o  and 
the  error  in  course  was  less  than  10°  in  83'- o  of  the  different  geometries. 

2.  With  a  small  bearing  error  (standard  deviation  =  0.5°),  only  33.4%  of  the  runs  had 
an  error  in  speed  of  less  than  10%  .  In  15.r'o  of  the  runs,  the  error  in  course  was 
less  than  10°. 

One  of  the  conclusions  of  this  study  is  that  the  Bearing  Extrapolation  Procedure,  even 
as  modified  here,  does  not  produce  accurate  estimates  of  Target  course  and  speed.  With 
Target  range,  however,  more  success  was  obtained.  Table  5  shows  the  percentage  of 
simulation  runs  that  produced  an  estimated  Target  range  at  the  end  of  the  maneuver 
within  5"o.  I0%s  and  20%>  of  the  true  range. 
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Table  5.     ACCURACY        IN        ESTIMATED        RANGE        WHEN        ONLY 
LAG-LEG-ESTIMATION  IS  USED.     (LEG  LENGTH  =   10  MIN). 


Own  Ship 
speed 
(knts) 

St.  Dev. 

bearing  er- 
ror 

Initial  dis- 
tance 

Percentage  of  times  the  absolute  error  is 
smaller  than  5%,  lO^.o,  and  20*^0  of  ac- 
tual for  range. 

5^0 

\o% 

20% 

15 

1.5° 

large 

21.6 

41.7 

72.1 

medium 

28.3 

54.6 

82.2 

short 

10.9 

19.7 

40.9 

1.0° 

large 

30.6 

57.1 

84.6 

medium 

35.4 

62.9 

86.6 

short 

10.7 

18.6 

40.2 

0.5° 

large 

49.9 

78.8 

93.7 

medium 

43.8 

71.1 

88.6 

short 

10.3 

18.9 

39.6 

0.0° 

large 

72.2 

88.9 

94.4 

medium 

55.6 

77.8 

88.9 

short 

11.1 

22.2 

38.9 

10 

1.5° 

large 

15.2 

30.6 

55.9 

medium 

20.2 

39.4 

70.7 

sliort 

16.4 

32.4 

51.3 

1.0° 

large 

21.0 

41.2 

70. 1 

medium 

26.2 

49.5 

76.6 

short 

16.1 

32.^ 

52.6 

0.5° 

large 

35.4 

62.1 

88.7 

medium 

35.4 

61.3 

79.5 

short 

15.3 

32.3 

52.5 

0.0° 

large 

72.2 

77.8 

94.4 

medmm 

38.9 

72.2 

77.8 

short 

16.7 

27.8 

50.0 

From  Table  5  it  can  be  seen  that  the  larger  the  Own  Ship  speed  ,  the  better  is 
the  estimate  of  range.  This  is  shown  in  Figure  9  which  is  the  plot  of  the  accuracy  in 
estimated  range  versus  standard  deviation  of  bearing  errors  when  the  Own  Ship  speed 
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is  10  and  15  knots.    This  accuracy  is  the  percentage  of  times  that  the  absolute  error  is 
within  20%  of  the  actual  for  range  at  the  end  of  the  maneuver. 
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Figure  9.        Comparison  of  Accuracy  for  Different  Own  Ship  Speeds  by  Using  only 
the  Lag-Leg-Estimation. 

Figure  10  shows  the  comparison  of  the  results  of  these  suggested  improvements  with  the 
results  obtained  in  Chapter  II  for  the  central  idea.  It  appears  that  for  short  initial  dis- 
tances the  accuracy  of  the  procedure  decreases,  but  for  large  and  medium  distances,  even 
with  large  bearing  error  (standard  deviation  =  1.5°),  the  accuracy  increases  consider- 
ably.  Note  that  the  accuracy  for  short  distances  was  not  adequate  in  either  case. 
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Figure  10.       Comparison  of  the  Accuracy  Between  Results  by  Using  only  the  Lag- 
Leg-Estimation  and  the  Results  by  Using  the  Central  Idea. 

3.     With  Previous  Minor  Changes  Plus  an  Increase  in  Turn  Duration 

In  addition  lo  the  two  minor  changes  made  before,  here  the  turn  duration  was 
increased  to  ten  minutes  and  the  simulation  experiments  were  repeated.  Table  6  shows 
the  results  for  this  case. 
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Table  6.  ACCURACY  IN  ESTIMATED 
LAG-LEG-ESTIMATION  IS  USED. 
Tl'RN  TIME  =   lOMIN.) 


RANGE        WHEN        ONLY 
(LEG   LENGTH    =    10  MIN. 


Own  Ship 
speed 
(knts) 

St.  Dev. 

bearing  er- 
ror 

Initial  dis- 
tance 

Percentage  of  times  the  absolute  error  is 
smaller  than  5%,  10%,  and  20%  of  ac- 
tual for  range. 

so/ 

10% 

2<)'''o 

15 

1.5° 

large 

25.5 

47.7 

78.3 

medium 

30.6 

54.3 

83.3 

short 

12.1 

21.7 

38.6 

1.0° 

larce 

34.8 

62.0 

87.8 

mediimi 

37.3 

59.9 

86.7 

short 

13.3 

22.9 

38.2 

0.5° 

large 

44.3 

73.5 

91.3 

medium 

40.7 

62.1 

87.8 

short 

14.4 

24.2 

38.6 

0.0° 

large 

66.7 

88.9 

94.4 

meditmi 

50.0 

61.1 

88.9 

short 

16.7 

22.2 

38.9 

10 

1.5° 

large 

18.3 

36.2 

65.1 

medium 

24.6 

46.2 

72.0 

short 

14.7 

28. 7 

50.3 

1.0° 

large 

25.9 

48.2 

78.4 

medium 

29.6 

53.9 

75.5 

short 

13.3 

28.1 

50.6 

0.5° 

large 

31.S 

59.0 

85.5 

mediuni 

32.2 

58. 6 

76.9 

short 

13,2 

26.6 

50.9 

0.0° 

large 

66.7 

77.8 

94.4 

medium 

66.7 

77.8 

short 

16.7 

-)-»  2 

44.4 

By  comparing  the  results  obtained  for  the  accuracy  in  range  in  Table  5  (Own 
ship  turn  =  6  mmutes)  and  Table  6  (Own  ship  turn  =  10  minutes),  it  is  appears  that 
small  improvements  in  precision  are  obtained  by  using  ten  minutes  for  the  turn. 
Figure   1 1  is  a  graphic  representation  of  this  fact. 
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Figure  11.       Comparison  of  Accuracy  for  Different  Own  Ship  Turn  Times  by  Using 
only  the  Lag- Leg- Estimation. 

a.     Effect  of  Target  Courses  Greater  than  I5(f  when  Lag-Leg  Estimation  and 
Minor  Changes  are  Used 

It  was  stated  in  Chapter  II  that  the  bearing  rate  for  the  lead-legs  does  not 
remain  constant  when  the  Target  course  is  greater  than  150°  (initial  angle  on  the  bow 
less  than  30°).  This  is  because  range  during  the  maneuver  decreases  to  the  point  where 
the  lead-leg  bearings  do  not  maintain  linearity.  To  determine  the  effect  on  the  esti- 
mation of  Target  range  that  this  behavior  produces,  simulations  were  conducted  as  be- 
fore except  that  results  were  obtained  for  12  separate  Target  courses  between  030°  and 
175°.  Previous  results  were  averaged  over  Target  course,  so  the  effect  of  a  particular 
Target  course  was  not  observed. 

The  results  of  one  simulation  experiment  are  depicted  in  Figure  12  and 
presented  in  Table  7.   The  parameters  used  to  obtain  these  results  were  the  following: 

1.  Initial  distance:  30,000  yds. 

2.  Target  couses:  030°,  050°,  070°,  090°,  110°,  130°,  150°,  155°,  160°,  165°,  170°,  and 

175°. 

3.  Target  speeds:  25,  15,  and  5  knots. 
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4.  Own  Ship  speed:  15  knots. 

5.  Standard  deviation  of  bearing  error:  1.5' 

6.  Leg  length:    10  minutes. 

7.  Turn  time:  10  minutes. 
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Figure   12.      Fraction  of  Times  Absolute  Error  is  within  20%  of  Actual  for  Range,  for 
Different  Target  Courses  and  Speeds. 
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Table  7.  TRACTION  OT  TIMES  ABSOLUTE  ERROR  IS  WITHIN  20%  OF 
ACTUAL  FOR  RANGE,  FOR  DIFFERENT  TARGET  COURSES  AND 
SPEEDS. 


Target 
speed 
(knts) 

Target  course  (°) 

030 

050 

070 

090 

110 

130 

150 

155 

160 

165 

170 

175 

25 

0.87 

0.83 

0.72 

0.98 

1.00 

1.00 

0.78 

0.62 

0.00 

0.00 

0.00 

0.00 

15 

0.86 

0.95 

0.95 

0.97 

0.99 

1.00 

0.94 

0.62 

0.38 

0.26 

0.13 

0.05 

5 

0.94 

0.96 

0.94 

0.97 

0.95 

0.93 

0.93 

0.94 

0.91 

0.97 

0.91 

0.95 

From  Figure  12  and  Table  7.  it  can  be  seen  that  for  a  Target  speed  of  5 
knots,  approximately  equal  accuracy  was  obtained  for  all  Target  courses.  But  for  Target 
speeds  of  15  and  25  knots,  much  less  accuracy  was  obtained  when  the  Target  course  was 
between  150°  and  175°.  For  Target  courses  less  than  150°,  there  was  less  observed 
change  in  the  estimation  of  Target  range. 

C.     USE  OF  AN  ARCTANGENT  MODEL  TO  EXTRAPOLATE  THE  LAG-LEG 
BEARINGS 

As  presented  so  far,  a  major  problem  with  the  Bearing  Extrapolation  procedure  is 
that  the  lag-leg  bearings  cannot  be  linearly  extrapolated  forward  or  backward  in  time 
while  maintaining  accuracy.  The  solution  to  this  problem  suggested  in  the  previous  se- 
veral sections  was  simply  not  to  use  extrapolated  bearings  from  the  lag-leg.  The  problem 
with  this  approach  is  that  tracking  information  is  lost  when  these  extrapolated  bearings 
are  not  used.  In  this  section,  an  attempt  will  be  made  to  recover  these  extrapolated 
bearings  by  using  an  arctangent,  rather  than  linear,  extrapolation  model. 

It  is  known  that  without  change  in  Own  Ship  course,  bearing  to  a  constant  course 
and  speed  moving  Target  is  of  the  form 


5,  =  5o  +  tan  '[ 


Vn 


R 


it-to)l 


(3.0) 


cpa 


where 
Bn   is  the  bearing  to  the  Target  at  CPA2, 


2  Closest  I'ornt  of  Approach 
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Vfj  is  the  relative  speed  between  Target  and  Own  Ship, 

R,^^  is  the  range  to  the  Target  at  CPA, 

^0   is  the  time  at  which  CPA  is  to  occur.    [Ref.  8:  pp.  2] 

Assuming  that  the  extrapolated  bearings  for  the  lag-leg  follow  this  arctangent 
model,  two  difierent  attempts  were  made  in  order  to  include  these  bearings  in  the  sol- 
ution of  the  problem. 

1.     Fitting  the  Arctangent  Model  by  Using  Least  Squares 

An  attempt  was  made  to  use  a  least  squares  technique  to  fmd  the  values  of  ^o. 
-r-^,  and  /()  in  (3.0)  which  produced  the  best  fit  to  the  observed  bearing  data.  Although 
this  approach  was  ultimately  unsuccessful,  there  may  be  another  least  squares  procedure 
that  will  prove  more  robust.  To  help  any  subsequent  researcher,  the  details  of  this  at- 
tempt are  presented  to  illustrate  what  has  already  been  examined. 
Let 


D 

^cpa 


Now  make  the  change  of  variable 

}>tan(/?,-5o), 


or 


Then  applying  least  squares,  it  is  necssar}'  to  minimize 


n 


A  necessarv  condition  is 


^  =  -2y ( };  -  a/  +  y.tQ)i  +  lYi}]  -a.i  +  a/o)'o  =  0,  (3.1) 


ca  i_^     '  "  ^^ 

:=1  i=\ 


and 
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n 

ess 


^  =  2^{ };  -  a/  +  a/o)a  =  0.  (3.2) 

°  i=\ 

Since  a  cannot  be  zero  unless  relative  speed  (Fq)  is  zero,  from  (3.2), 


i=\ 


^{Y,-at  +  at,)  =  0.  (3.3) 


Dividing  (3.3)  by  n  and  solving  for  a. 


K=  -  ^      .  (3.4) 

r-/n 


Now  substituting  (3.3)  in  (3.1)  gives 


n 

YiY^-a)  1^  +  ^1X1  =  ^.  (3.5) 

tl^  ^  7^ 


Substituting  (3.4)  in  (3.5)  and  solving  for  t^  yields 


n 

n 
;=1 

/     J         ' 

^  =  -^^ ^ •  (3-6) 


Because  the  determination  of  }',  for  all  /,  requires  knowledge  of  the  value  of  the 
bearing  at  CPA  (^0),  an  iterative  procedure  was  attempted  where  the  first  value  for  Bq 
was  guessed,  then  with  equations  (3.5)  and  (3.6),  a  and  ^0  were  computed.  Then  by  using 


;=1 
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a  new  value  for  B^  was  obtained.  This  iteration  was  continued  in  an  attempt  to  converge 
on  the  true  B^,.  Unfortunately  for  even  small  errors  in  bearings,  (standard  deviation  = 
0.5°).  the  procedure  would  not  reliably  converge. 

2.     Estimation  of  the  Parameters  for  the  Arctangent  Model  by  Solving  a  System  of 
Linear  Equations 

It  is  known  that  with  three  difierent  bearings  to  a  Target  that  is  moving  with  a 
constant  course  and  speed,  the  unknown  parameters  of  the  equation 

5,  =  i5o  +  tan~'[a(/-/o)] 

can  be  determined.  [Ref  7:  pp.  2] 

Now  let  5,.  B,.  and  Z?3  be  the  bearings  obtained  at  times  ^,,  z^.  ^rid  /,  respec- 
tively. Then  it  can  be  written 

B,  =  B^  +  \^n-\cf.{t,-iQ)-\, 

B2  =  5o  +  tan-'[a(/2-ro)]. 

53  =  Bo  +  tair'[o;(/3-/o)]. 

An  attempt  will  be  made  to  solve  for  B^,  a.,  and  /"p.    Substracting, 

B.-  B^=  tan"'[a(/2  -  /q)]  -  tan"'[a(/,  -  /q)], 

B^-  B,=  tan"^[(x(/3  -  /q)]  -  tan~'[a(/i  -  /q)], 

B^-B2  =  tan"'[a(/3  -  /q)]  -  tan"' [7.(^2  -  /q)]- 

Now  taking  the  tangent  of  both  sides  and  knowing  the  trigonometric  identity  for  the 
tangent  of  the  difference  of  two  angles,  the  following  can  be  obtained: 


tanf5:,  -  /?,) 


<h  -  h~) 


1  +  aV:  -  ^oK^i  -  ^) 


tan(i53  -  i?,j  = ^ 


1  +  a  (/3  -  /u)(^i  -  ^0) 
tan(53  -/?-,)  = ^ — 
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From  these  follows  that 


1  taniB.-B,) 

ih  -  h)(A  -  h)  tan(52  -  ^i)  =  -f  (/2  -  h) '-\ '-  ,  (3.7) 

a 

1  tan(B.-B,) 


1                     tan{B.  -  B-/) 
(/3-/o)(/2-ro)tan(53-52)  =  -F(^3-^2) H -•  (^-9) 

Dividing  (3.7)  by  (3.8).  (3.7)  by  (3.9),  and  (3.8)  by  (3.9),  allows  the  solution  for  t^,  a,  and 

[.  tan(/?3  -  B.)  tan(Z?2  -  B^)Uh  -  t^)i[j  -  /,)  -  (^3  -  /,)(/2  -  ^,)]  +  /,(/i  -  ^2)^ 
^0  =  - 


(/2  -  ^3)  tan(52  -  5,)[(/3  -  ^1)  tan(^3  -  Bj)  -  (h  -  (2)  tan(^3  -  5,)] 
tan(Z?2  -  B^)  tan(Z?3  -  Z?,)  +  /3(/2  -  i^f  tan(53  -  B2)  tan(Z?3  -  5,) 
-  (^2  -  /,)  tan(53  -  BMh  -  h)  tan(52  -  B,)  -  (t^  -  t,)  iMB^  -  5,)] 

{[.  -  I2)  taniB.  -  B{)  tan(i?3  -  5,) 
^~  {^2-  h){h  -  ^1)  tan(/?2  -  5,)  +  (^2  -  ^i)('o  -  h)  tan(53  -  B^)  ' 

i?o=5i-tan"'[a(/i-;o)]. 

A  subroutine  to  compute  these  values  was  written  and  the  simulation  of  the 
problem  was  done  with  the  following  conditions: 

1.  Speed  of  Own  Ship:    5.  10,  and  15  knts. 

2.  Speed  of  the  Target:    5,  15,  and  25  knts. 

3.  Initial  range:    10.000,  30.000,  and  60,000  yds. 

4.  Initial  angle  in  the  bow  for  Own  Ship:    90°. 

5.  Legs  length:    10  ininutes. 

6.  Time  bet\\een  legs:    6  minutes. 

7.  Standard  deviation  bearing  error:   0°,  0.5°,  1.0°,  and  1.5°. 
S.   Only  lead-lag  maneuver  is  considered. 

9.  Bearings  for  the  first  leg  were  faired  and  extrapolated  linearly. 

10.  Bearings  in  the  second  leg  were  faired  by  using  a  quadratic  model. 

11.  'Ihe  extrapolation  of  the  bearings  for  the  second  leg  was  done  by  using  the  for- 
mulas obtained  above  and  by  using  faired  bearings  from  that  leg.    This  was  done 
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because  the  use  of  raw  bearings  produced  too  much  error  in  the  estimation  of  the 
parameters  of  the  Target. 

Results  are  presented  in  Table  8.    It  is  emphasized  that  these  results  are  only 
for  lead-lag  maneuvers. 


Table  8.  ACCURACY  IN  ESTIMATED  RANGE  BY  USING  THE 
ARCTANGENT  EXTRAPOLATION  MODEL  FOR  THE  LAG-LEG 
(LEAD-LAG  MANELrV'ER). 


Own  Ship 
speed 
(knts) 

St.  Dev. 
bearing  er- 
ror 

Initial  dis- 
tance 

Percentage  of  times  the  absolute  error  is 
smaller  than  5%,  10%,  and  20%  of  ac- 
tual for  range. 

S% 

lo-^o 

20% 

15 

1.5° 

large 

20.4 

40.4 

70.6 

medium 

26.0 

50.9 

77.8 

short 

11.5 

25.2 

46.9 

1.0° 

large 

2^.6 

54.8 

83.5 

medium 

31.1 

59.2 

83.8 

short 

11.0 

25.2 

47.1 

0.5° 

large 

45.9 

75.5 

93.4 

medium 

37.S 

66.8 

87.2 

short 

17.2 

31.8 

52.1 

0.0° 

large 

66.7 

88.9 

94.4 

medium 

55.6 

66.7 

88.9 

short 

11.1 

27.8 

50.0 

10 

1.5° 

large 

16,4 

30.9 

58.6 

mediuni 

23.0 

43.8 

77.1 

short 

19.6 

34.4 

59.4 

1.0° 

large 

22.9 

43.0 

72.5 

medium 

29.  S 

53.9 

83.9 

short 

19.5 

34.4 

60.8 

0.5° 

large 

34.9 

60.3 

87.2 

medium 

32.9 

59.1 

78.2 

short 

l.v2 

2^M 

56.  1 

0.0° 

large 

7-1  ■) 

77.8 

Q4.4 

medium 

38.9 

72.2 

77.8 

short 

11.1 

27.8 

55.6 
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Comparing  these  results  with  those  on  Table  3  on  page  13,  it  is  clear  that  this 
procedure  performs  better  than  the  original  (unimproved)  Bearing  Extrapolation  Proce- 
dure. However,  the  results  are  about  the  same  as  those  in  Table  8,  indicating  that  sim- 
ply not  using  the  lag-leg  extrapolated  bearings  is  just  about  as  efiective  as  trv'ing  to 
extrapolate  them  wdth  the  arctangent  model. 
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IV.     CONCLUSIONS  AND  RECOMMENDATIONS 

1.  CONCLUSIONS 

1.  The  proposed  procedure  does  not  allow  estimation,  with  adequate  accuracy,  of  the 
Target  track  parameters  when  linear  bearing  extrapolation  is  used  for  both  legs 
(central  idea).  Linear  bearing  extrapolation  is  only  valid  for  lead-legs  where  low 
and  constant  bearing  rate  is  obtained.  For  lag-legs  the  bearing  rate  is  generally 
higher  and  more  variable. 

2.  Adequate  accuracy  for  course  and  speed  is  not  obtained  using  any  of  the  examined 
procedures  unless  bearing  error  is  zero  (ideal  conditions). 

3.  Adequate  accuracy  in  estimated  range  for  large  and  medium  initial  distances  can 
be  obtained  if: 

a.  Quadratic  faired  bearings  from  the  lag-leg  are  used  in  combination  with  linear 
extrapolated  bearings  from  the  lead  leg. 

b.  Lead-lag  maneuver  is  performed. 

c.  Lead  angle  close  to  90°  is  used. 

d.  Leg  length  is  larger  than  seven  minutes  and  less  than  twelve  minutes. 

e.  Own  ship  turn  time  is  larger  than  six  minutes  and  less  than  ten  minutes. 

f    Angle  on  the  bow  is  not  smaller  than  30°.  If  it  is,  then  Target  speed  is  not 
greater  than  5  knots. 

4.  Adequate  accuracy  in  range  for  short  initial  distances  is  obtained  only  when  the 
Target  is  not  closing  the  Own  Ship. 

2.  RECOMMENDATIONS 

1.  If  this  procedure  is  implemented,  use  only  extrapolated  bearings  from  the  lead  leg 
(low  and  constant  bearing  rate)  and  faired  bearings  for  the  lag-leg  (high  and  vari- 
able bearing  rate). 

2.  Further  research  should  be  done  to  obtain  a  model  that  allows  use  of  the  extrapo- 
lated bearines  obtained  from  the  las-les. 
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APPENDIX  A.     DISCUSSION  OF  SELECTED  SOLUTION  METHODS 

A.  SPIESS  RANGING 

Spiess  theory  establishes  that  the  problem  of  obtaining  satisfactory  solutions  for  a 
bearings-only  approach  against  targets  moving  with  a  constant  course  and  speed,  can 
be  solved  by  using  four  bearings  and  at  least  one  change  of  course  or  speed  by  Own 
Ship.  With  this  in  mind,  several  graphical  methods  capable  of  giving  quick  reliable  sol- 
utions were  developed.   [Ref  3] 

From  the  Passive  Ranging  Manual,  Volume  III,  the  basic  graphical  method  for 
Spiess  ranging  is  explained  as  follows: 

Given  three  bearings  to  the  Target,  observed  at  times  r,,  /j  »  ^rid  /j  together 
with  a  fourth  time,  /j.  the  locus,  L.  of  all  possible  Target  positions  at  time,  t^,  for  all 
possible  Target  tracks  which  satisfy  the  three  bearings  at  /,,  t^,  and  i^,  and  which 
maintain  a  constant  course  and  speed,  is  a  straight  line.  The  actual  position  of  the 
Target  at  time,  z^.  is  determined  by  the  intersection  of  L  with  the  observed  bearing 
at  time,  u,  (provided  SSK  changes  course  and  speed). 

This  locus.  L,  can  be  determined  by  picking  out  two  arbitran."  Target  tracks 
which  satisfy  the  three-bearing  conditions,  and  plotting  the  position  of  the  Target 
on  these  tracks  at  time  i^  Then  L  is  the  line  through  these  two  plotted  positions. 
This  is  the  basis  of  the  Spiess  four-bearing  T.VIA.  The  bearing  lines  at  t^.  I2  ,  and  t^ 
are  used  to  find  the  locus  of  Target  positions  at  i^.  The  intersection  of  this  locus 
with  the  bearing  hne  at  t^  determines  one  point  on  the  track.  Another  point  is  de- 
termined by  using  the  bearing  lines  at  /,.  ^3  ,  and  t^  to  find  the  locus  of  Target  posi- 
tions at  /i-  The  intersection  of  this  locus  with  the  bearing  line  at  /,  determines 
another  point  on  the  track.  The  track  is  then  the  line  joining  the  two  constructed 
Target  positions. 

Because  bearings  are  not  precise  enough  to  apply  this  method,  faired  bearings  are 
usually  used. 

B.  EKELUND  RANGING 

The  basic  Ekelund  solution  for  a  constant  course  and  speed  moving  Target  is  based 
on  a  two  leg  maneuver  by  Own  Ship.  During  these  legs.  Own  Ship  records  actual  Target 
bearings.  The  rate  of  change  of  the  bearings  is  computed  for  each  leg  by  applying  linear 
least  squares  techniques.  This  method  assumes  that  during  the  turn  the  distance  that 
Own  Ship  moves  and  the  change  of  angle  in  the  bow  can  be  disregarded;  i.e.,  assumes 
an  instantaneous  turn.  Given  this  assumption,  the  range  to  the  Target  can  be  computed. 
Using  Figure  13.  derivation  of  range  is  as  follows: 
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leg  one 


target  track 


Figure  13.      Geometry  for  EkeJund  Ranging  Method. 


The  velocity  of  the  Target  across  the  line  of  sound  for  Own  Ship's  first  leg  is 


{A.\) 


and  for  leg  2 


V,n=V^smd-,  =  R^^f-  =  R,Ba 


{A.2) 


Own  ship's  velocity  across  the  line  of  sound  for  leg  1  can  be  expressed  as 


cB, 


V,,,=  V,,sm4>,  =  R,-jf-=R,B,, 


{A.3) 


and  for  leg  2 


J^.o2=^o2Sin(/>2  =  /?2^=^242 


(^.4) 


where 
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B„  =  Bearing  rale  of  Target  in  leg  i. 
/?!  =  Range  to  the  Target  at  the  turn. 
R2  =  Range  to  the  Target  after  the  turn. 
F„,  =  Velocity  of  Own  Ship  in  leg  i. 
V„  =  Velocity  of  the  Target  in  leg  i. 

Assuming  Target  velocity  remains  constant  during  Own  Ship  maneuvering;  i.e., 
V,^  =  r',2  =  V,  and  that  the  angle  in  the  bow  also  remains  approximately  constant;  i.e., 
0^  =  0^  =  6,  and  that  the  range  to  the  Target  does  not  appreciably  changes;  i.e., 
R^  =  R2  =  R„  then  equations  (A.l),  (A. 2),  (A. 3),  and  (A.4)  can  be  written  as 

V,sme^=R,B,^,  {A. 5) 

lysine  =  R,B^,  {A. 6) 

and  by  subtracting  (A. 7)  from  (A. 5)  and  (A. 8)  from  (A. 6) 

V,  sm  e  -  V,,  sm  0,  =  R^B,.  -  B,,)  =  R,B„ 

\\  sm  e  -  V,2  sm  02  =  7^,(^2  "  ^2)  =  ^e^2- 
Subtracting  these  last  two  equations  yields 

r;,  sin(/)i  -  V^2  sin  02 


Re  = 


Bi-B, 


where,  B,  =  bearing  rate  of  Target  relative  to  Own  Ship  in  leg  i. 
The  problems  in  accuracy  with  this  method  are: 

1.  It  assumes  that  Own  Ship  can  perform  the  two  legs  maneuver  with  an  instantane- 
ous turn  [Ref  9;  pp.  1-3]. 

2.  The  Ekelund  Range  Equation  is  dependent  on  the  ratio  of  the  bearing  rates  devel- 
oped on  the  two  legs  of  the  ranging  maneuver  which  are  small  and  not  ver>"  precise 
values  [Ref  9:  pp.  1-3]. 
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3.  In  practice,  measuring  the  bearing  rate  involves  a  bearing  smoothing  process  which 
takes  place  over  a  portion  of  each  leg  and  results  in  bearing  rates  at  times  signif- 
icantly different  from  the  time  of  the  turn  [Ref  9:  pp.  3-2]. 

C.     CHURN  METHOD 

The  basic  concept  of  CHURN  TMA  is  similar  to  that  of  the  strip  plot  and  includes 
in  its  solution  the  Spiess  TMA  for  the  particular  case  where  there  are  only  four  bearing 
lines  computed  [Ref  9:  pp.  7-1,  7-3].  Figure  14,  depicts  the  regular  CHURN  TMA 
method. 


0siimated 
target  track 


Bi  =  bearing 
at  time  I 


Figure   14.      Geonietr>  for  CHURN  Method. 

Own  Ship  performs  a  two  leg  maneuver  during  which  it  measures  bearings  to  a 
Target  that  is  moving  with  constant  course  and  speed.  The  CHURN  method  takes  those 
bearings  and  fits  a  Target  track  to  them  using  the  perpendicular  distance  between  a  given 
bearing  line  and  the  corresponding  Target  position  on  the  fitted  track  as  a  measure  of 
the  goodness  of  fit.  Applying  the  concept  of  least  squares,  the  best  Target  track  is  that 
which  then  minimizes  the  sum  of  squares  of  these  distances. [Ref  9:  pp.  7-2  ] 
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APPENDIX  B.     VARIABLES  OF  THE  SIMULATION  PROGRAM 

1.  A.B.C.D.E.F.  Entries  of  the  A''^  matrix  for  quadratic  regression  in  the  first  leg. 
These  values  were  computed  by  using  APL  [Ref.  10].  Then  they  were  introduced 
as  initial  data. 

2.  A0.A1.A2.   CoefTicients  for  the  quadratic  model  in  the  first  leg. 

3.  AOl.Al  1.A22.   Coefilcients  for  the  quadratic  model  in  the  second  leg. 

4.  B22.  Temporary  variable  used  to  compute  the  estimated  position  of  the  target  on 
the  estimated  target  track  at  the  begining  of  the  maneuver. 

5.  BEERR.   Bearing  error. 

6.  BEEXTRdX).    Vector  of  extrapolated  bearings. 

7.  BEE.  Temporar}'  variable  used  to  compute  the  estimated  position  of  the  target  on 
the  estimated  target  track  at  the  end  of  the  maneuver. 

8.  BE(IX).   Vector  of  bearings  to  the  target. 

9.  BEIIX.   Y  intercept  for  the  linear  model  of  bearings  versus  time  in  the  first  leg. 

10.  BEl(lX).   Temporar\'  vector  of  bearings  to  the  target. 

11.  BE21X.   Y  intercept  for  the  linear  model  of  bearings  versus  time  in  the  second  leg. 

12.  BMEAX.    VIean  of  the  bearing  to  the  target.    Same  variable  is  used  for  both  legs. 

13.  BO.    Bearing  to  the  target  at  CPA. 

14.  CC.    Counter  used  in  the  computation  of  the  precision  of  the  procedure. 

15.  CIXCO.    Times  that  absolute  value  of  the  range  error  is  less  than  5%. 

16.  ClXCOl.   Times  that  absolute  value  of  the  speed  error  is  less  than  5%. 

17.  CIXC02.    Times  that  absolute  value  of  the  course  error  is  less  than  5°. 
IS.  C0EL(2).   Vector  of  own  ship  courses  for  the  first  leg. 

19.  CQSL.    Own  ship  course  for  the  second  leg. 

20.  CJ\   Percentage  of  times  that  error  in  estimated  target  course  is  less  than  5°. 

21.  CIO.    Percentage  of  times  that  error  in  estimated  target  course  is  less  than  10°. 

22.  C20.    Percentage  of  times  that  error  in  estimated  target  course  is  less  than  20°. 

23.  DEGPl.    Pi  expressed  m  degrees. 

24.  DIEZ.    Times  that  absolute  value  of  the  range  error  is  less  than  \0%. 

25.  DIEZl.   Times  that  absolute  value  of  the  speed  error  is  less  than  10%. 

26.  D1EZ2.    Times  that  absolute  value  of  the  course  error  is  less  than  10°. 

27.  DSEED.    Seed  for  the  pseudorandom  generation  of  normal  bearing  errors. 
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28.  ECTLR.   Estimated  target  course  obtained  by  linear  regression  on  the  target  fixes. 

29.  ER.ANTL.    Estimated  range  to  the  target  at  the  end  of  the  maneuver. 

30.  ESLOPE.   Temporary  variable  corresponding  to  the  target  course. 

31.  ESPTL.    Estimated  speed  of  the  target. 

32.  EXT(IX).    Vector  of  estimated  target  positions  in  the  X  axis. 

33.  EYINT.   Y  intercept  of  the  estimated  target  course. 

34.  E^"T(TN).    Vector  of  estimated  target  positions  in  the  Y  axis. 

35.  FBE(IN).   Vector  of  faired  target  bearings. 

36.  IX.   Time  to  complete  the  maneuver  in  1/3  of  minutes. 

37.  IXFO(IN).  Flag  used  to  signal  when  Estimated  X  position  of  the  target  is  equal 
to  the  X  position  of  the  own  ship. 

38.  IX(3).   Vector  of  initial  target  distances  from  own  ship  in  the  X  axis. 

39.  n"(3).    Vector  of  initial  target  distances  from  own  ship  in  the  Y  axis. 

40.  LLL.  Counter  to  determine  the  amount  of  fixes  used  in  the  regression  to  determine 
target  track. 

41.  PER.  Percentage  error  in  range  computed  from  the  estimated  and  the  simulated 
range  to  the  target  at  the  end  of  the  maneuver. 

42.  PERI.    Percentage  error  in  speed  computed  from  the  estimated  and  the  simulated 

target  speed. 

43.  PER2.  Amount  of  error  in  course  computed  from  the  estimated  and  the  simulated 
target  course. 

44.  PL    Pi  in  radians. 

45.  R.S.T.L.V.W.  Entries  of  the  /i"'  matrix  for  quadratic  regression  in  the  second  leg. 
These  values  were  computed  by  using  APL  [Ref.  lOj.  Then  they  introduced  as  ini- 
tial data. 

46.  R.ANGEF.    Simulated  range  to  the  target  at  the  end  of  the  maneuver. 

47.  RCBFLl.  Rate  of  change  of  bearing  during  the  first  leg.  This  value  is  computed 
assuming  the  linear  model  for  the  bearings  versus  time. 

48.  RCBSLl.  Rate  of  change  of  bearing  during  the  Second  leg.  This  value  is  computed 
assuming  the  linear  model  for  the  bearings  versus  time. 

49.  R5.   Percentage  of  times  that  error  in  estimated  range  to  the  target  is  less  than  5%. 

50.  RIO.    Percentage  of  times  that  error  in  estimated  ranee  to  the  tareet  is  less  than 

l()"o. 

51.  R20.    Percentage  of  times  that  error  in  estimated  range  to  the  target  is  less  than 

20^0. 

52.  SCT(6).    Vector  of  simulated  target  courses. 

53.  SD.    Standard  deviation  of  beanne  errors. 
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54.  SLOPE(I\).   Vector  of  faired  bearings  to  the  target  expressed  as  slopes  in  the  X-Y 
plane. 

55.  SLOPEX(IX).    Vector  of  extrapolated  bearings  to  the  target  expressed  as  slopes 
in  X-Y  plane. 

56.  SMALL.  Constant  that  defines  equality  between  simulated  Own  ship  position  and 

estimated  target  position. 

57.  SMALLB.    Constant  used  to  determine  when  the  faired  or  extrapolated  target 
bearing  correspond  to  a  value  of  90°  or  270°. 

58.  S0(3).    Vector  of  simulated  own  ship  speeds. 

59.  SST(3).   Vector  of  simulated  target  speeds. 

60.  SL.VIB.   Summation  of  the  bearing  values. 

61.  SLMB2.   Summation  of  the  square  of  the  bearing  values. 

62.  SL'.VIT.    Summation  of  time. 

63.  SUMTB.    Summation  of  the  products  of  bearing  and  time. 

64.  SLMT2.    Summation  of  the  squares  of  the  time. 

65.  SU\1T2B.    Summation  of  the  products  of  square  of  time  and  bearing. 

66.  SUMX.    Summation  of  the  estimated  target  positions  in  the  X  axis. 

67.  SUMXY.    Summation  of  the  product  of  estimated  target  positions  in  the  X  axis 
and  the  estimated  target  position  in  the  Y  axis. 

68.  SUMX 2.   Summation  of  the  square  of  the  estimated  target  positions  in  the  X  axis. 

69.  SLMY.    Summation  of  the  estimated  target  positions  in  the  ^'  axis. 

70.  SL"M^'2.   Summation  of  the  square  of  the  estimated  target  positions  in  the  Y  axis. 

71.  SXT(I\).   Vector  of  simulated  target  positions  in  the  X  axis. 

72.  S'^"T(1N).   Vector  of  simulated  target  positions  in  the  Y  axis. 

73.  S5^  Percentage  of  times  that  error  in  estimated  target  speed  is  less  than  5%. 

74.  SIO.    Percentage  of  times  that  error  in  estimated  target  speed  is  less  than  10%. 

75.  S20.    Percentage  of  times  that  error  in  estimated   target  speed  is  less  than  20%. 

76.  TEL(2).    Vector  of  leg  lengths. 

77.  TEMP.    Temporan."  variable. 

78.  TMEAN.    Mean  of  the  time. 

79.  TO.    1  ime  when  CPA  occurs. 

80.  TWO.    2  expressed  as  double  precision  constant. 

81.  VARERB.    Variance  of  bearing  error. 

82.  VEINTE.    Times  that  absolute  value  of  the  range  error  is  less  than  20%. 
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83.  VEINTl.   Times  that  absolute  value  of  the  speed  error  is  less  than  20%. 

84.  VEINT2.   Times  that  absolute  value  of  the  course  error  is  less  than  20°. 

85.  VD.    Ratio  between  relative  speed  of  target  and  own  ship  and  range  at  CPA. 

86.  XQ(IN).   Vector  of  own  ship  positions  in  the  X  axis. 

87.  XOEXTR(IX).    Vector  of  extrapolated  own  ship  positions  in  the  X  axis. 

88.  XVI EAN.    Mean  of  target  positions  (fixes)  in  the  X  axis. 

89.  XEE.   Estimated  target  position  at  the  end  of  the  maneuver  in  the  X  axis. 

90.  X22.    Estimated  target  position  at  the  begining  of  the  maneuver  in  the  X  axis. 

91.  "^'FE.    Estimated  target  position  at  the  end  of  the  maneuver  in  the  Y  axis. 

92.  YMEAN.    Mean  of  target  positions  (fixes)  in  the  Y  axis. 

93.  YO(IN).    Vector  of  own  ship  positions  in  the  Y  axis. 

94.  YQEXTR(IN).   Vector  of  extrapolated  own  ship  positions. 

95.  Y22.    Estimated  target  position  at  the  beggining  of  the  maneuver  in  the  Y  axis. 
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APPENDIX  C.     SIMULATION  PROGRAM  CODE 

PROGRAM  SIMTMA 

ji^  ^<.  ^^  y*  .J.  .j^  k*.  ^<«  .*<>•..*.  .t<  .*..*.  ^'.  .ji.  .'^  ^Y  .jT.  ^*«  ^.  y.  ^^  y^  ^^  J- ^^  .^^ 

"  THIS  PROGRAM  SIMULATES  A  PASSIVE  BEARINGS-ONLY  TARGET  MOTION  ANALYSIS''^ 
•-  PROCEDURE  WHICH  IS  PERFORMED  BY  A  SUBMARINE  (OWN  SHIP)  AGAINST  A  * 
^-  TARGET  WHICH  IS  MOVING  WITH  CONSTANT  COURSE  AND  SPEED.  * 

"  THE  OWN  SHIP  PERFORMS  A  TWO  LEG  MANEUVER.  EACH  OF  THEM  WITH  CONSTANT  * 
''-  COURSE  AND  SPEED.  * 

i;  ie    it  Vf  Vr         i:    i(  Vr  ic  it 

it  it 

*  WRITTEN  BY  CMR.  CARRERO  CUBEROS,  BERNABE  * 

^'  VENEZUELAN  NAVY  * 

of.  ^  oU  .j^  ^^  .J.  ji^  ^Y  -*'  ^(  V'  "*'  <^'  ^'  "^c  "*'  ^c  i'  'if  '^'  Vf  V'  ^'c  y*  ^V  "*'  V'  Vc  Vr  Vr  Vc  V^  Vc  Vc  Vr  "ic  Vc  Vr  V'  V'  Vc  Vc  V'  Vf  Vc  Vc  Vr  V'  "'c  ^'  V?  ^c  V'  V'  Vc  Vc  V?  Vc  V'  Vc  Vf  ^p  V'  Vc  Vc  ic  V'  V'  Vc  Vc  V? 

PARAMETERdN  =78) 

INTEGER  TEL(2),  INFO(IN),  LLL,  CC ,  CINCO,  DIEZ,  VEINTE, 
&CINC01,  DIEZl,  VEINTl,  CINC02,  DIEZ2,  VEINT2 

DOUBLE  PRECISION  IX(3),  IY(3),  DSEED,  C0FL(2),  COSL,  SCT(6),  SD, 
&SST(3),  VARERB,  BEERR,  EXT(IN),  EYT(IN),  S0(3),  SXT(IN),  SYT(IN), 
&XO(IN),  YO(IN),  BE(IN),  SLOPE(IN),  BEEXTR(IN),  XOEXTR(IN), 
&YOEXTR(IN),  SLOPEX(IN),  PI,  SMALL,  DEGPI ,  TWO,  TEMP,  SUMX,  SUMX2 , 
&SUMY2,  SUMY,  SUMXY,  XMEAN,  YMEAN,  ESLOPE,  SMALLB ,  EYINT,  ECTLR, 
&B22,  BFF,  X22,  Y22,  XFF,  YFF,  ESPETL,  ERANTL,  RANGEF,  SUMT,  SUMB, 
&SUMT2,  SUMB2,  SUMTB ,  TMEAN,  BMEAN,  RCBFLl,  RCBSLl,  BEl(IN),  BEIIN, 
&BE2IN,  FBE(IN),  AO ,  Al,  A2 ,  SUMT2B ,  A,  B,  C,  D,  E,  F,  AOl,  All, 
&A21,  R,  S,  T,  U,  V,  W,  BO,  TO,  VD,  PER,  PERI,  PER2 ,  R5 ,  RIO,  R20, 
&.SS,    SIO,  S20,  C5,  CIO,  C20 

DATA  PI/3.  1415926535898D0/,  SMALL/0.  IDO/ ,  DEGPI/180. DO/ ,  TWO/2. DO/, 
&SCT/153.  DO  ,  158.  DO  ,  163.  DO  ,  167.  DO  ,  171.  DO  ,  175.  DO/  ,  TEL/39  ,  30/  , 
&SO/30000. DO, 20000. DO, 10000. DO/ ,  SST/50000. DO, 30000.  DO, 10000.  DO/ , 
&COFL/088.D0,272.D0/,  SMALLB/. OOOOIDO/ ,  A/ 118. 4785317D0/ , 
&B/-3. 77185762D0/,  C/. 02947719689D0/ ,  D/. 1205858692D0/ , 
&E/-. 0009459915779D0/,  F/. 000007448752582D0/ ,  R/. 3438423645D0/ , 
&S/-. 0450738916D0/,  T/. 001231527094D0/ ,  U/. 007603190052D0/ , 
&V/-. 00023091133D0/,  W/. 000007448752582D0/ , 
&IY/60000. DO, 30000. DO, 10000. DO/,  IX/0. D0,0. D0,0. DO/ 

PRINT-'S' ENTER  VARIANCE  OF  BEARING  ERROR  (SQUARE  DEGREES)  AND' 

PRINT-"-, 'INITIAL  SEED  (7  DIGIT  INTEGER)  TO  BEGIN  THE  SIMULATION' 

READ--'-,  VARERB,  DSEED 

•^'   TRANSFORMATION  OF  DATA  THAT  INVOLVES  DEGREES  TO  RADIANS 
SD  =  SQRT( VARERB) 
SD  =  SD^'-PI/DEGPI 
DO  50  I  =  1,2 

COFL(I)  =  COFL(I)--'^PI/DEGPI 
50    CONTINUE 

DO  100  I  =  1,6 

SCT(I)  =  SCT(I)-->PI/DEGPI 
100    CONTINUE 
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•>''   SELECTION  OF  TWO  DIFFERENT  LENGTH  LEG 
DO  2000  I  =  1,2 

^^  SELECTION  FROM  THREE  DIFFERENT  INITIAL  DISTANCES  BETWEEN  TARGET  AND 

*  OWN  SHIP 

DO  1999  MMMM  =1,3 

*  SELECTION  FROM  THREE  DIFFERENT  SPEEDS  FOR  OWN  SHIP 

DO  1100  L  =  1,3 

''-  SELECTION  FROM  LEAD-LAG  AND  LAG-LEAD  MANEUVER  AND  INITIALIZATION 
DO  950  M  =  1,1 
CC  =  0 
CINCO  =  0 
DIEZ  =  0 
VEINTE  =  0 
CINCO 1  =  0 
DIEZl  =  0 
VEINTl  =  0 
CINC02  =  0 
DIEZ2  =  0 
VEINT2  =  0 

*  SELECTION  FROM  SIX  DIFFERENT  TARGET  COURSES 

DO  900  J  =  1,6 

•^-   SELECTION  FROM  THREE  TARGET  SPEEDS 
DO  850  K  =  1,3 

*  REPETITION  OF  THE  SIMULATION  100  TIMES  WITH  DIFFERENT  SEED 

DO  793  JIJ  =  1,100 

^^  INITIALIZATION  AND  SIMULATION  OF  TARGET  POSITION  EVERY  20  SECONDS 

DO  150  II  =  1,2'-TEL(I) 

SXT(II)  =  O.DO 

SYT(II)  =  0.  DO 

XO(II)  =  O.DO 

YO(II)  =  O.DO 

BE(II)  =  O.DO 

BEl(II)  =  O.DO 

FBE(II)  =  O.DO 

SLOPE(II)  =  O.DO 

SLOPEX(II)  =  O.DO 

XOEXTR(II)  =  O.DO 

YOEXTR(II)  =  O.DO 

BEEXTR(II)  =  O.DO 

INFO(II)  =  O.DO 

TEMP  =  (DBLE(II))/180.D0 

SXT(II)  =  IX(MMMM)  +  COS(SCT(J)-PI/TWO)^-SST(K)''--TEMP 

SYT(II)  =  lY(MMMM)  -  SIN(  SCT(  J) -PI/T\\'0)''^SST(K) -TEMP 
150  CONTINUE 

-'•  INITIALIZATION  OF  REGRESSION  TO  COMPUTE  RATE  OF  CHANGE  OF  BEARING  IN 
•'■  FIRST  LEG. 

SUMT  =  O.DO 

SUMT2  =  0.  DO 
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SUMB2  =  O.DO 

SUMB  =  O.DO 

SUMTB  =  O.DO 

SUMT2B  =  O.DO 

DO  750  N  =  l,TEL(I)-9 

*  SIMULATION  OF  POSITION  OF  OWN  SHIP  EVERY  20  SECONDS  IN  THE  FIRST  LEG. 

TEMP  =  (DBLE(N))/180.D0 

XO(N)  =  COS(COFL(M)  -  PI/TWO)*SO(L)*TEMP 

YO(N)    =  SIN(PI/T\v'0    -   COFL(M))'VSO(L)'>TEMP 

*  SIMULATION  OF  THE  BEARINGS  TO  THE  TARGET  IN  THE  FIRST  LEG. 

CALL  NORRN(DSEED,BEERR) 
IF(ABS(SYT(N)  -  YO(N))  .  GE.  SMALL)THEN 
IF(ABS(XO(N)  -  SXT(N))  . GE.  SMALL)THEN 
IF(YO(N)  .LT.  SYT(N))THEN 
IF(XO(N)  .  LT.  SXT(N))THEN 

BE(N)  =  ATAN((SXT(N)  -  XO(N) )/(SYT(N)  -  YO(N))) 
ELSE 
TEMP  =  ATAN((SXT(N)  -  XO(N) )/(SYT(N)  -  YO(N))) 
BE(N)  =  TWO''--PI  +  TEMP 
END  IF 
ELSE 
IF(XO(N)  .LT.  SXT(N))THEN 
TEMP  =  ATAN((SXT(N)  -  XO(N) )/(SYT(N)  -  YO(N))) 
BE(N)  =  PI  +  TEMP 
ELSE 
TEMP  =  ATAN((SXT(N)  -  XO(N) )/(SYT(N)  -  YO(N))) 
BE(N)  =  PI  +  TEMP 
END  IF 
END  IF 

BE(N)  =  BE(N)  +  BEERR*SD 

IF(BE(N)  .GE.  TWO''^PI)BE(N)  =  BE(N)  -  TWO''^PI 
ELSE 
IF(SYT(N)  .GT.  YO(N))THEN 
BE(N)  =  O.DO 
ELSE 

BE(N)  =  PI 
END  IF 
INFO(N)  =  99 
END  IF 
ELSE 
IF(SXT(N)  .GT.  XO(N))THEN 

BE(N)  =  PI/TWO 
ELSE 

BE(N)  =  3.D0'>PI/TW0 
END  IF 
END  IF 

'■>-  COMPUTATION  OF  COEFFICIENTS  FOR  QUADRATIC  AND  LINEAR  LEAST  SQUARES 

*  FITS  FOR  THE  BEARINGS  IN  THE  FIRST  LEG. 

SUMT  =  SUMT  +  DBLE(N) 

SUMT2  =  SUMT2  +  (DBLE(N)  )''">2 

IF(BE(N)  .GT.  3.D0--'--PI/TW0)THEN 

BEl(N)  =  BE(N)  -  TU^0--'^PI 
ELSE 
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BEl(N)  =  BE(N) 
END  IF 

SUMS  =  SUMB  +  BEl(N) 
SUMB2  =  SUMB  2  +  BEl(N)*''f2 
SUMTB  =  SUMTB  +  (DBLE(N) )*BE1(N) 
SUMT2B  =  SUMT2B  +  (DBLE(N)''">2)''^BE1(N) 
CONTINUE 

TMEAN  =  SUMT/DBLE(TEL(I)-9) 
BMEAN  =  SUMB/DBLE(TEL(I)-9) 
RCBFLl  =  (SUMTB  -  SUMT''-BMEAN)/(SUMT2  - 
BE  UN  =  BMEAN  -  RCBFL1''^TMEAN 
AOl  =  R''^SUMB  +  S-'^SUMTB  +  T''^SUMT2B 
All  =  S'^SUMB  +  U-'^SUMTB  +  V'^SUMT2B 
A21  =  T'^SUMB  +  V"  SUMTB  +  W--'^SUMT2B 


SUMT^'^TMEAN) 


*  COMPUTATION  OF  THE  OWN  SHIP  POSITION  AT  THE  FICTITIOUS  TURN 

TEMP  =  DBLE(TEL(I))/180.D0 

XO(TEL(I))  =  COS(COFL(M)  -  PI/TWO)*SO( L)*TEMP 

YO(TEL(I))  =  SIN(PI/TWO  -  COFL(M)  )'VSO(L)''^TEMP 

*  COMPUTATION  OF  FAIR  BEARINGS  DURING  THE  FIRST  LEG  USING  THE  RESULTS 

*  FROM  LINEAR  OR  QUADRATIC  FIT.  THIS  DEPEND  IF  THE  MANEUVER  IS  LEAD-LAG 

*  OR  LAG -LEAD. 

DO  755  N  =  l,TEL(I)-9 
IF  (M  .EQ.  1)THEN 

FBE(N)  =  RCBFLl'^^DBLECN)  +  BEIIN 
ELSE 

FBE(N)  =  AOl  +  All''fDBLE(N)  +  A21''KDBLE(N)*''^2) 
END  IF 

IF(FBE(N)  .LT.  O.DO)FBE(N)  =  TWO^'^PI  +  FBE(N) 
IF((ABS(FBE(N)  -  PI/T^'O)  .  LT.  SMALLB)  .OR. 
&  (ABS(FBE(N)  -  3.D0'>PI/TVv'0)  .  LT.  SMALLB ))THEN 

SLOPE(N)  =  O.DO 
ELSE 

SLOPE(N)  =  l.DO/TAN(FBE(N)) 
END  IF 
755  CONTINUE 

'^  IN  CASE  OF  LAG-LEAD  MANEUVER  (M=2),  CALL  SUBROUTINE  TANGENT  TO  COMPUTE 
'^  THE  COEFFICIENTS  THAT  RELATES  THE  BEARINGS  TO  THE  ARCTANGENT  MODEL.  IN 
■>''   CASE  OF  LEAD-LAG  MANEUVER  USE  THE  LINEAR  MODEL  ALREADY  COMPUTED.  WITH 
■>''   THIS  INFORMATION,  COMPUTE  THE  EXTRAPOLATED  BEARINGS  FROM  THE  FIRST  LEG 

IF(M  .EQ.  2)CALL  TNGENT(FBE , IN,M,BO,TO , VD) 
DO  760  N  =  TEL(I)+10,2*TEL(I) 
IF(M  .EQ.  1)THEN 

BEEXTR(N)  =  RCBFLl^^DBLE(N)  +  BEIIN 
ELSE 

BEEXTR(N)  =  BO  +  ATAN( VD*(DBLE(N)  -  TO)) 
END  IF 

IF(BEEXTR(N)  .  LT.  0.  D0)BEEXTR(N)=TVv'O^>PI+BEEXTR(N) 
IF((ABS(BEEXTR(N)  -  PI/TWO)  . LT.  SMALLB)  .OR. 
&  (ABS(BEEXTR(N)  -  3.  DO'-'-PI/TWO)  .  LT.  SMALLB))THEN 

SLOPEX(N)  =  O.DO 
ELSE 

SLOPEX(N)  =  l.DO/TAN(BEEXTR(N)) 
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END  IF 
760  CONTINUE 

*  SIMULATION  OF  EXTRAPOLATED  POSITIONS  OF  OWN  SHIP  OUT  OF  THE  FIRST  LEG. 

DO  765  N  =  TEL(I)+10,  2*TEL(I) 
TEMP  =  (DBLE(N))/180.D0 

XOEXTR(N)  =  COS(COFL(M)  -  PI/TWO)*SO(L)*TEMP 
YOEXTR(N)  =  SIN( PI/TWO  -  COFL(M) )*SO(L)*TEMP 

765  CONTINUE 

*  SIMULATION  OF  FAIR  BEARING  TO  THE  TARGET  AT  THE  TURN,  BY  USING  THE 

*  LINEAR  MODEL. 

FBE(TEL(I))  =  RCBFL1*DBLE(TEL(I))  +  BEIIN 
IF(FBE(TEL(I))  . LT.  0. DO)FBE(TEL( I) )  =  TWO*PI  + 

&  FBE(TEL(I)) 

IF(FBE(TEL(I))  .  GE.  TWO'>PI)FBE(TEL(  I) )  =FBE(TEL(I))  - 

&  TWO*PI 

*  DETERMINATION  OF  OWN  SHIP  COURSE  FOR  THE  SECOND  LEG  ACCORDING  WITH  THE 

*  LINE  OF  SIGHT  TO  THE  TARGET  AT  THE  TURN. 

IF  (M  .EQ.  1)THEN 

COSL  =  FBE(TEL(I))  -   PI'H  • 48888889D0) 
ELSE 

COSL  =  FBE(TEL(I))  +  PI'>( .  48888889D0) 
END  IF 

IF(COSL  .LT.  O.DO)COSL  =  TWO^^PI  +  COSL 
IF(COSL  .GT.  TWO''>-PI)COSL  =  COSL  -  TWO''-PI 

'"^  INITIALIZATION  OF  REGRESSION  TO  COMPUTE  RATE  OF  CHANGE  OF  BEARING  IN 
''^  SECOND  LEG. 

SUMT  =  O.DO 

SUMT2  =  0.  DO 

SUMB2  =  O.DO 

SUMB  =  O.DO 

SUMTB  =  O.DO 

SUMT2B  =  0.  DO 

DO  778  N  =  TEL(I)+10,2*TEL(I) 

*  SIMULATION  OF  POSITION  OF  OWN  SHIP  EACH  20  SECONDS  IN  THE  SECOND  LEG. 

TEMP  =  (DBLE(N)  -  DBLE(TEL( I ) ) )/ 180.  DO 

XO(N)  =  XO(TEL(I))  +  SO(L)'''TEMP-'"^COS(COSL-PI/TWO) 

YO(N)  =  YO(TEL(I))  -  SO(L)''^TEMP*SIN(COSL-PI/TWO) 

''^  SIMULATION  OF  THE  BEARINGS  TO  THE  TARGET  IN  THE  SECOND  LEG. 

CALL  NORRN(DSEED,BEERR) 
IF(ABS(SYT(N)  -  YO(N))  .  GE.  SMALL)THEN 
IF(ABS(SXT(N)  -  XO(N))  .  GE.  SMALL)THEN 
IF(Y0(N1  . LT.  SYT(N))THEN 
IF(X0(N)  .LT.  SXT(N))THEN 

BE(N)  =  ATAN((SXT(N)  -  XO(N) )/( SYT(N)  -  YO(N))) 
ELSE 
TEMP  =  ATAN((SXT(N)  -  XO(N) )/( SYT(N)  -  YO(N))) 
BE(N)  =  TWO-'>PI  +  TEMP 
END  IF 
ELSE 
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IF(XO(N)  .LT.  SXT(N))THEN 
TEMP  =  ATAN((SXT(N)  -  XO(N) )/(SYT(N) 
BE(N)  =  PI  +  TEMP 
ELSE 
TEMP  =  ATAN((SXT(N) 


YO(N))) 


XO(N))/(SYT(N)  -  YO(N))) 


+  BEERR'-^SD 
O.DO)BE(N) 


BE(N)  =  PI  +  TEMP 
END  IF 
END  IF 

BE(N)  =  BE(N) 
IF(BE(N)  .LT. 
IF(BE(N)  .GE. 
ELSE 
IF(SYT(N)  .GT. 
BE(N)  =  O.DO 
ELSE 

BE(N)  =  PI 
END  IF 
INFO(N)  =  99 
END  IF 
ELSE 
IF(SXT(N) 

BE(N) 
ELSE 

BE(N) 
END  IF 
END  IF 


=  TWO'^PI  + 


TW'0'>PI)BE(N)  =  BE(N) 
YO(N))THEN 


BE(N) 
■  TWO'-^PI 


GT.  XO(N))THEN 
=  PI/TWO 

=  3.D0'"^PI/TW0 


*  COMPUTATION 

*  BEARINGS  IN 


778 


OF  COEFICIENTS  FOR  LINEAR  AND  QUADRATIC  FITS  FOR  THE 
THE  SECOND  LEG 

SUMT  =  SUMT  +  DBLE(N) 

SUMT2  =  SUMT2  +  (DBLE(N)  )*''^2 

IF(M  .EQ.  1)THEN 

IF(BE(N)  .GT.  .65D0'>PI)THEN 

BEl(N)  =  BE(N)  -  TWO'''PI 
ELSE 

BEl(N)  =  BE(N) 
END  IF 
ELSE 

BEl(N)  =  BE(N) 
END  IF 

SUMB  =  SUMB  +  BEl(N) 
SUMB2  =  SUMB  2  +  BE1(N)^">2 
SUMTB  =  SUMTB  +  (DBLE(N)  )''^BE1(N) 
SUMT2B  =  SUMT2B  +  (DBLE(N)-'-''^2)^'^BEl(N) 
CONTINUE 

TMEAN  =  SUMT/DBLE(TEL(I)-9) 
BMEAN  =  SUMB/DBLE(TEL(I)-9) 

RCBSLl  =  (SUMTB  -  SUMT''-BMEAN)/(  SUMT2  -  SUMT^'TMEAN) 
BE2IN  =  BMEAN  -  RCBSL1''^TMEAN 
AO  =  A^'^SUMB  +  B^'^SUMTB  +  C^-SUMT2B 
Al  =  B^'^SUMB  +  D^'^SUMTB  +  E''^SUMT2B 
A2  =  C^'-SUMB  +  E^--SUMTB  +  F''SUMT2B 


COMPUTATION  OF  FAIR  BEARINGS  DURING  THE  SECOND  LEG  USING  LINEAR 
OR  QUADRATIC  FIT.  THIS  DEPEND  ON  THE  TYPE  OF  MANEUVER. 
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TEL(I)+10,2*TEL(I) 

DTHEN 
=  AO  +  A1*DBLE(N)  + 


779 


DO  779  N  = 
IF(M  .EQ 

FBE(N) 
ELSE 

FBE(N)  =  RCBSL1*DBLE(N)  +  BE2IN 
END  IF 

IF(FBE(N)  .LT.  0 
IF((ABS(FBE(N)  - 

(ABS(FBE(N)  -  3.D0*PI/TW0) 

SLOPE(N)  =  O.DO 
ELSE 

SLOPE(N)  = 
END  IF 
CONTINUE 


A2*(DBLE(N)**2) 


DO)FBE(N)  =  TWO*PI  +  FBE(N) 
PI /TWO)  .LT.  SMALLB)  .OR. 

LT.  SMALLB)) THEN 


l.DO/TAN(FBE(N)) 


*  IN  CASE  OF  LEAD-LAG  MANEUVER  (M=l),  CALL  SUBROUTINE  TANGENT  TO  COMPUTE 

*  THE  COEFFICIENTS  THAT  RELATES  THE  BEARINGS  TO  THE  ARCTANGENT  MODEL.  IN 

*  CASE  OF  LAG-LEAD  MANEUVER  USE  THE  LINEAR  MODEL  ALREADY  COMPUTED.  WITH 

*  THIS  INFORMATION,  COMPUTE  THE  EXTRAPOLATED  BEARINGS  FROM  THE  SECOND  LEG 

IF(M  .EQ.  DCALL  TNGENT(FBE,  IN,M,BO,TO,VD) 
DO  781  N  =  l,TEL(I)-9 
IF(M  .EQ.  DTHEN 

BEEXTR(N)  =  BO  +  ATAN(VD'V(DBLE(N)  -  TO)) 
ELSE 

BEEXTR(N)  =  RCBSL1*DBLE(N)  +  BE2IN 
END  IF 

*  WRITE(1,'(F9. 6,1X,I3)')BEEXTR(N),N 

IF(BEEXTR(N)  .  LT.  0.  D0)BEEXTR(N)=TWO'VPI+BEEXTR(N) 
IF((ABS(BEEXTR(N)  -  PI/TWO)  .  LT.  SMALLB)  .OR. 

&  (ABS(BEEXTR(N)  -  3. DO*PI/TWO)  . LT.  SMALLB) )THEN 

SLOPEX(N)  =  O.DO 
ELSE 

SLOPEX(N)  =  l.DO/TAN(BEEXTR(N)) 
END  IF 
781  CONTINUE 


*   SIMULATION  OF  EXTRAPOLATED  OWN  SHIP  POSITIONS  OUT  OF  THE  SECOND  LEG. 

DO  782  N  =  l,TEL(I)-9 
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TEMP  =  (DBLE(TEL(I))  -  DBLE(N) )/180.  DO 

XOEXTR(N)  =  XO(TEL(I))  -  COS(COSL  -  PI/TWO)* 

& 

SO(L)'^TEMP 

YOEXTR(N)  =  YO(TEL(I))  +  SIN(COSL  -  PI/T\s'0)'V 

6c 

SO(L)''^TEMP 

CONTINUE 

''•-  INITIALIZATION  OF  VARIABLES  USED  IN  THE  LINEAR  REGRESSION  TO  COMPUTE 
■>''   COURSE  OF  THE  TARGET. 

SUMX  =  O.DO 

SUMX2  =  0.  DO 

SUMY2  =  O.DO 

SUMY  =  O.DO 

SUMXY  =  0.  DO 

LLL  =  0 

'•-  COMPUTATION  OF  EXPECTED  POSITION  (X,Y)  OF  THE  TARGET. 
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DO  783  N  =  1,2*TEL(I) 

*  ELIMINATION  OF  18  CENTRAL  ESTIMATED  POSITIONS  OF  THE  TARGET  TO  ALLOW 

*  THE  MANEUVER  OF  THE  OWN  SHIP  IN  ORDER  TO  CHANGE  COURSE  FROM  FIRST  LEG 

*  TO  SECOND  LEG  (SIX  MINUTES) 

IF((N  .GT.  TEL(I)-9)  .AND.  (N  .  LE.  (TEL(I)+9))) 
&  GO  TO  783 

IF(INFO(N)  .EQ.  99)THEN 

EXT(N)  =  XO(N) 

EYT(N)  =  SLOPEX(N)*(XO(N)  -  XOEXTR(N))  +  YOEXTR(N) 
ELSE 

EXT(N)  =  XOEXTR(N) 

EYT(N)  =  SLOPE(N)''f(XOEXTR(N)  -  XO(N))  +  YO(N) 
END  IF 

*  COMPUTATION  OF  VALUES  OF  THE  VARIABLES  NEEDED  TO  DO  LINEAR  REGRESSION 

*  ON  THE  TARGET  FIXES  TO  DETERMINE  ESTIMATED  COURSE 

LLL  =  LLL  +  1 

SUMX  =  SUMX  +  EXT(N) 

SUMX2  =  SUMX2  +  (EXT(N)  )''^''^2 

SUMY  =  SUMY  +  EYT(N) 

SUMY2  =  SUMY2  +  (EYT(N)  )'^''^2 

SUMXY  =  SUMXY  +  EXT(N)'VEYT(N) 
783  CONTINUE 

XMEAN  =  SUMX/ LLL 
YMEAN  =  SUMY /LLL 

*  ESTIMATED  SLOPE  (ESLOPE)  AND  Y  INTERCEPT  (EYINT)  OF  TARGET  TRACK 

ESLOPE  =  (SUMXY  -  SUMX->YMEAN)/( SUMX2  -  SUMX-^XMEAN) 
EYINT  =  YMEAN  -  ESLOPE-'^XMEAN 

'•f  ESTIMATION  OF  COURSE,  SPEED  AND  RANGE  OF  THE  TARGET 

ECTLR  =  (PI/TU'O  -  ATAN( ESLOPE ))''>-DEGPI/PI 

B22  =  EYT(l)  -  EYINT  -  ESLOPE''EXT(  1) 

BFF  =  EYT(2-->TEL(I))  -  EYINT  -  ESLOPE''-EXT(  2''^TEL(  I) ) 

Y22  =  EYT(l)  -  B22/(l  +  ESLOPE--"'-- 2) 

YFF  =  EYT(2"-TEL(I))  -  BFF/(1  +  ESL0PE''"'^2) 

X22  =  (Y22  -  EYINT)/ESL0PE 

XFF  =  (YFF  -  EYINT) /ESLOPE 

ESPETL  =  SQRT((XFF  -  X22)-'">2+( YFF  -  Y22)''"'^2)/ 

&  DBLE(2->TEL(I)-1) 

ERANTL  =  SQRT((XFF  -  X0(  2*TEL(  I)  )  )''^''^2  + 

&  (YFF  -  Y0(2''^TEL(I)))''"'^2) 

RANGEF  =  SQRT((SXT(2''^TEL(I))  -  X0(  2'>TEL(  I) )  )'"f^^2  + 

6c  ( SYT(  2'-TEL(  I ) )  -  Y0(  2'^TEL(  I) )  )*"2) 

*  INITIAL  COMPUTATIONS  TO  DETERMINE  THE  PRECISION  OF  THE  PROCEDURE 

CC  =CC  +  1 

'^  PRECISION  FOR  ESTIMATED  TARGET  RANGE 

PER  =  ((ERANTL  -  RANGEF)*100. DO)/RANGEF 
IF(ABS(PER)  .LE.  20.  DO)THEN 
VEINTE  =  VEINTE  +  1 

IF(ABS(PER)  .LE.  10.D0)THEN 
DIEZ  =  DIEZ  +  1 
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IF(ABS(PER)  .LE.  5.D0)THEN 

CINCO  =  CINCO  +  1 
END  IF 
END  IF 
ELSE 
END  IF 

*  PRECISION  FOR  ESTIMATED  TARGET  SPEED 

PERI  =  ((ESPETL^nSODO  -  SST(K)  )*100.  DO)/SST(K) 
IF(ABS(PER1)  .LE.  20.D0)THEN 
VEINTl  =  VEINTl  +  1 

IF(ABS(PER1)  .LE.  10.D0)THEN 
DIEZl  =  DIEZl  +  1 

IF(ABS(PER1)  .LE.  5.  DO)  THEN 

CINCOl  =  CINCOl  +  1 
END  IF 
END  IF 
ELSE 
END  IF 

*  PRECISION  FOR  ESTIMATED  TARGET  COURSE 

PER2  =  ECTLR  -  SCT( J)''^DEGPI/PI 
IF(ABS(PER2)  .LE.  20.  DO)THEN 
VEINT2  =  VEINT2  +  1 

IF(ABS(PER2)  .LE.  10.  DO)THEN 
DIEZ2  =  DIEZ2  +  1 

IF(ABS(PER2)  .LE.  5.D0)THEN 

CINC02  =  CINC02  +  1 
END  IF 
END  IF 
ELSE 
END  IF 

*  OUTPUT  INDICATING  SIMULATED  AND  ESTIMATED  TARGET  TRACK  PARAMETERS 

WRITE(1,'(I3,2X,F9.  1,1X,F12.  1 , 2X, 2F5.  1 , IX, 2F11.  1) 

&  ')M,RANGEF,ERANTL,SCT(J)'''DEGPI/PI, ECTLR, SST(K), 

S.  ESPETL'nSO.  DO 

793  CONTINUE 

795  CONTINUE 

850  CONTINUE 

900  CONTINUE 

*  FINAL  COMPUTATIONS  TO  DETERMINE  THE  PRECISION  OF  THE  PROCEDURE. 

R5  =  (DBLE(CINCO)/DBLE(CC))-'nOO 
RIO  =  (DBLE(DIEZ)/DBLE(CC))''aOO 
R20  =  (DBLE(VEINTE)/DBLE(CC))''-100 
S5  =  (DBLE(CINC01)/DBLE(CC))^nOO 
SIO  =  (DBLE(DIEZl)/DBLE(CC))'nOO 
S20  =  (DBLE(  VEINTl) /DBLE(CC))--aOO 
C5  =  (DBLE(CINC02)/DBLE(CC))--nOO 
CIO  =  (DBLE(DIEZ2)/DBLE(CC))"100 
C20  =  (DBLE(VEINT2)/DBLE(CC))*100 

''-  OUTPUT  THAT  SHOWS  THE  PRECISION  OF  THE  PROCEDURE 

WRITE(1,'(I1,1X,F7. 1,1X,F3. 1,1X,F8. 1,2X,9F5. 1)') 
&  M,  SO(L),  (SD/PI)''-DEGPI,  lY(MMMM),  R5 ,  RIO,  R20, 
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&  S5,  SIO,  S20,  C5,  CIO,  C20 

950       CONTINUE 
1100     CONTINUE 

1999  CONTINUE 

2000  CONTINUE 
STOP 

END 

*  THIS  SUBROUTINE  GENERATES  ONE  NORMAL  (0,1)  NUMBER  BY  USING  THE  BOX  * 

*  MULLER  METHOD.  THIS  NUMBER  IS  USED  BY  THE  MAIN  PROGRAM  TO  CREATE   * 

*  THE  ERRORS  IN  THE  BEARINGS  WHICH  ARE  ASSUMED  NORMAL  (0,VARERB)     * 

SUBROUTINE  NORRN(DSEED,BEERR) 

DOUBLE  PRECISION  U(2),  D31M1,  DSEED,  TWO,  PI,  BEERR 

DATA  D31M1  /2147483647. DO/ ,  TWO/2. DO/,  PI/3.  141592654D0/ 

*  LINEAR  CONGRUENTIAL  METHOD  TO  OBTAIN  A  UNIF0RM(0,1) 

DO  5  I  =1,2 

DSEED  =  DMOD(  16807.  D0'-DSEED,D31M1) 
U(I)  =  DSEED/D31M1 
5     CONTINUE 

*  BOX  AND  MULLER  METHOD  TO  OBTAIN  A  NORMAL(0,1) 

BEERR  =  (SQRT(  -  TWO*LOG(U(  1) ) )  )*COS(TWO'''PI*U(  2)  ) 

RETURN 
END 

■>''   SUBROUTINE  TO  COMPUTE  THE  PARAMETERS  OF  THE  ARCTANGENT  MODEL  (BO, 

*  TO,  AND  VD)  BY  USING  FAIRED  BEARINGS  FROM  THE  QUADRATIC  MODEL.       •>'' 

*  NOTE:  THIS  CASE  IS  USEFUL  ONLY  FOR  TEN  MINUTE  LEGS.  'V 

SUBROUTINE  TNGENT( FBE , IN , M , BO , TO , VD) 
PARAMETER(PI  =  3.  141592654D0) 

DOUBLE  PRECISION  BO,  TO,  VD,  XY,  YZ ,  ZA,  Tl ,  T2 ,  T3 ,  FBE(IN), 
6cDE,  EF,  FG,  GH,  SUMTO ,  SUMVD,  SUMBO 
INTEGER  AB,  BC,  CD,  M,  IN,  COUNT,  JK,  KJ 
SUMTO  =  O.DO 
SUMVD  =  0.  DO 
SUMBO  =  O.DO 
COUNT  =  0 
IF(M  .EQ.  1)THEN 

JK  =  49 

KJ  =  78 
ELSE 

JK  =  1 

KJ  =  30 
END  IF 
DO  6  III  =  JK,(KJ-14) 

AB  =  III 

BC  =111+7 

CD  =  III  +  14 

Tl  =  DBLE(AB) 

T2  =  DBLE(BC) 

T3  =  DBLE(CD) 

DE  =  T2  -  T3 

EF  =  T3  -  Tl 
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FG  =  T3  -  T2 

GH  =  T2  -  Tl 

XY  =  TAN(FBE(BC)  -  FBE(AB)) 

YZ  =  TAN(FBE(CD)  -  FBE(AB)) 

ZA  =  TAN(FBE(CD)  -  FBE(BC)) 

TO  =  (XY''^YZ*Tl''^(FG''f''^2)  -  ZA^^XY*T2*(EF**2)  +  ZA*YZ*T3'V(GH**2) )/ 
X       (DE'-XY''^(ZA-''^EF  -  YZ''^FG)  -  GH'>ZA'V(XY*EF  -  YZ''^GH)) 

VD  =  (XY''^YZ'>DE)/(XY''^EF'^(T2  -  TO)  +  YZ''^GH'>(TO  -  T3)) 

BO  =  FBE(AB)  -  ATAN(VD'V(T1  -  TO)) 

IF(BO  .GT.  2.D0*PI)B0  =  BO  -  2.  DO^PI 

IF(BO  .GT.  3.D0-'^PI/2.D0)B0  =  BO  -  2.  DO''-PI 

SUMTO  =  SUMTO  +  TO 

SUMVD  =  SUMVD  +  VD 

SUMBO  =  SUMBO  +  BO 

COUNT  =  COUNT  +  1 
CONTINUE 

TO  =  SUMTO/DBLE( COUNT) 
VD  =  SUMVD /DBLE( COUNT) 
BO  =  SUMBO/DBLE( COUNT) 
RETURN 
END 
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